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^ ABSTRACT 

Aims. We aim to constrain the properties and evolutionary status of early and mid-spectral type supergiants (from 04 to 07.5). These 
posses the highest mass-loss rates among the O stars, and exhibit conspicuous wind profiles. 

Methods. Using the non-LTE wind code CMFGEN we simultaneously analyzed the FUV-UV and optical spectral range to determine 
the photospheric properties and wind parameters. We derived effective temperatures, luminosities, surface gravities, surface abun- 

I dances, mass-loss rates, wind terminal velocities, and clumping filling factors. 

Results. The supergiants define a very clear evolutionary sequence, in terms of ages and masses, from younger and more massive 
stars to older stars with lower initial masses. 04 supergiants cluster around the 3 Myr isochrone and are more massive than 60 M a , 
while the 05 to 07.5 stars have masses in the range 50 - 40 M e and are 4±0.3 Myr old. The surface chemical composition is typical 
of evolved O supergiants (nitrogen-rich, carbon- and oxygen-poor). While the observed ranges of carbon and nitrogen mass-fractions 
are compatible with those expected from evolutionary models for the measured stellar masses, the N/C ratios as a function of age 
^3 are inconsistent with the theoretical predictions for the four earliest (04 spectral type) stars of the sample. We question the efficiency 

Oh °f rotational mixing as a function of age for these stars and suggest that another mechanism may be needed to explain the observed 

abundance patterns. Mass-loss rates derived with clumped-models range within a factor of three of the theoretical mass-loss rates. 
The corresponding volume-filling factors associated with small-scale clumping are 0.05±0.02. Clumping is found to start close to the 

> photosphere for all but three stars, two of which are fast rotators. 

d Key words. Stars: winds - Stars: atmospheres - Stars: early-type - Stars: fundamental parameters 



> 1. Introduction 

in 

£ — , Many fields of modern astrophysics are directly or indirectly re- 
lated to the physics of massive stars. Locally, massive stars are 
C<*) responsible for the ionization and expansion of their surround- 
1^ ing HII regions, and are believed to be responsible for triggering 
>) star formation (e.g. Zavagno et al. 2007; Martins et al. 2010). 
^vq More generally, they dominate the chemical and dynamical evo- 
y— \ lution of the global interstellar medium of their host galaxies. 
J* Starburst events, either in the local or the distant universe, are 
. dominated by hot massive stars, as revealed by the composite 
spectrum of distant, star-forming galaxies, which exhibit numer- 
S_i ous features typical of O and B stars (e.g. Steidel et al. 1996). 
Fast rotating massive stars are likely progenitors of some of the 
most energetic events in the Universe, such as core-collapse su- 
pernovae and gamma-ray-bursts (GRBs) (e.g. Hjorth et al. 2003; 
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Woosley & Bloom 2006) and population III stars are believed 
to be very massive objects with extreme effective temperatures 
(Bromm et al. 1999; Nakamura & Umemura 1999). 

In the 1970s it was recognized that stellar evolution in the 
upper Hertzsprung-Russell (H-R) diagram is, to a large extent, 
governed by mass loss. Through the course of their evolution, 
massive stars loose a significant fraction of their mass through 
their stellar winds (see e.g., the review by Meynet & Maeder 
2007), and through giant eruptions (e.g. Smith & Owocki 2006). 
Although the line-driven wind theory has been successful in ex- 
plaining the global behavior of stellar winds, the actual wind 
parameters remain subject to uncertainties, especially regard- 
ing the wind terminal velocity, mass-loss rate and the degree 
of dumpiness (see e.g., the workshop held in 2007 in Potsdam, 
Hamann et al. 2008). The most striking uncertainty regarding 
the mass-loss rates is related to the intrinsically unstable nature 
of the line -driving mechanism (Lucy & Solomon 1970; Owocki 
et al. 1988), which leads to strong reverse shocks (e.g. Owocki 
& Puis 1999) that separate fast low-density wind material from 
overdense clumps. Since many prominent mass-loss diagnostics 
are sensitive to the square of the density, they overestimate the 
mass-loss rate in the presence of clumping. The uncertainty in 
the true value of the mass-loss rate translates in a corresponding 
uncertainty in the evolutionary tracks of massive stars. 
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Recent studies showed that winds are clumped, although 
there is no consensus concerning the amount, nature (optically 
thin versus optically thick) and stratification of clumping. Initial 
studies (Crowther et al. 2002; Hillier et al. 2003; Bouret et al. 
2003, 2005; Fullerton et al. 2006) found a reduction in mass- 
loss rates by a factor of three to ten when using the optically 
thin assumption for clumping. Within the same framework, the 
radial distribution of the clumping factor was constrained by 
Puis et al. (2006) and Najarro et al. (2011), based on a simul- 
taneous modeling of Ha infrared, millimeter and radio observa- 
tions. These studies concluded that clumping is three to six times 
stronger in the lower wind, where Ha forms, compared with 
the outer wind, where the radio continuum originates. These 
findings marginally agree with the theoretical predictions by 
Runacres & Owocki (2002, 2005). The aforementioned stud- 
ies assumed a void-interclump medium, although Zsargo et al. 
(2008) demonstrated that the treatment of a non-void interclump 
medium is required to explain the global spectroscopic proper- 
ties of massive stars. 

The optically thin approximation was subsequently chal- 
lenged by Oskinova et al. (2007), who claimed that significantly 
lower reduction (factor of three only) is found when using opti- 
cally thick clumping. When it is pronounced, this macroclump- 
ing (or porosity) affects the emerging spectrum for a given mass- 
loss by reducing the effective opacity of the wind, because op- 
tically thick clumps hide material. Another factor that might 
affect the mass-loss rate is the structure of winds in velocity 
space (hereafter called velocity porosity). This arises because 
the dense absorbing clumps along a given line of sight, only oc- 
cur at discrete velocities, and leads to a reduction in the strength 
of P Cygni absorption components (Owocki 2008). This velocity 
porosity has been implemented in 2D/3D stochastic wind mod- 
els by Sundqvist et al. (2010, 201 1), who concluded that it could 
increase mass-loss rates derived from UV resonance lines by up 
to an order of magnitude. 

Major progress has also been achieved in the last few years 
in stellar rotation theory, but more progress is needed. It has 
been known for a long time that the physical properties of stars 
are distorted by rotation, as predicted from the von Zeippel 
theorem and indicated from recent interferometric studies (e.g. 
Domiciano de Souza et al. 2003). Consequently, rotation is also 
expected to cause a latitudinal dependence of the line force and 
hence mass loss (Maeder & Meynet 2000). 

From the view point of stellar evolution, rotation causes an 
increase in the mass-loss rate, a change in evolutionary tracks in 
the H-R diagram, a lowering of the effective gravity, an extension 
of the main-sequence phase, the transport of angular momentum 
and the mixing of CNO-cycle processed material up to the stellar 
surface (Maeder & Meynet 2000; Heger & Langer 2000). The 
so-called helium discrepancy first identified by Herrero et al. 
(1992) in many galactic O-type stars and especially in fast ro- 
tators is naturally explained by evolution models with rotation 
— higher helium (and nitrogen) abundances are predicted for 
higher masses and higher rotation velocities. Another conse- 
quence of rotation is that two stars at the same evolutionary 
phase and at the same location in the H-R diagram can have 
different initial masses. 

Several other observational facts have been successfully ac- 
counted for using evolutionary models with rotation (Maeder & 
Meynet 2000). Models with rotation predict the correct blue-to- 
red supergiant ratio as observed in the Small Magellanic Cloud 
(Maeder & Meynet 2001), explain the observed population of 
Wolf-Rayet stars as a function of metallicity (Meynet & Maeder 
2003, 2005), and provide a consistent framework for the evo- 



lution of progenitors of soft-long GRBs (Yoon & Langer 2005; 
Woosley & Bloom 2006; Meynet & Maeder 2007). 

Despite these success, some observations remain unex- 
plained. Hunter et al. (2007, 2008, 2009), who investigated the 
link between surface abundance patterns and rotation in the 
Galaxy and in the Magellanic Clouds, concluded that up to 40% 
of their sample presented nitrogen surface abundances that could 
not be explained by rotational mixing. These studies strongly 
suggest that the global chemical evolution of massive stars, as 
probed from their surface abundances, is an intricate combina- 
tion of several physical parameters, including mass, metallicity, 
age (evolutionary status) and rotation (Maeder et al. 2009). 

Some of the spectral lines used to determine mass-loss rates 
and clumping factors originate from CNO ions. Because surface 
CNO abundance patterns are modified by rotational mixing dur- 
ing stellar evolution, it is very important to have reliable esti- 
mates of these patterns in order to fully characterize the wind 
properties of O-type stars. 

For this study, we have analyzed a sample consisting of eight 
Galactic supergiants with spectral types ranging from 04 If to 
07.5 If (Sota et al. 2011). We derive their properties, and ad- 
dress questions about as their evolutionary status and how their 
surface abundances relate to this. We also investigate whether 
a consistent picture for mass loss and clumping may be drawn 
from different diagnostics such as P v AA1 118-1 128, O v /U371, 
N iv /11718, He n /14686 and Ha. By using different diagnostics 
we can test the influence of the ionization structure, the adopted 
abundances, and to a limited extent, the clumping formulation 
on wind diagnostics. 

In the following section, we present our sample selection 
and observational material. The modeling tools are presented in 
Sect. 3, while diagnostics of the different stellar parameters are 
presented in Sect. 4. Results of the analysis for the individual 
objects are presented in Sect. 5 and a general discussion of the 
evolutionary status, chemical properties and wind properties is 
given in Sect. 6, before our general conclusion in Sect 7. 

2. Target selection and observations 

2.1. Stellar sample 

Table 1 presents the sample selected for this study and funda- 
mental observational data. The selected objects are considered 
to be representative of early-Of supergiants, whose surface prop- 
erties are likely to show processed material from stellar evolu- 
tion and whose strong winds are expected to exhibit conspicu- 
ous signatures of clumping. The spectral classification scheme 
of these stars was designed by Walborn (1971a,b, 1973) and was 
recently updated in Sota et al. (2011). The scheme is meant to 
sort stars with different temperatures, luminosities, and abun- 
dances. CNO-cycled material must play a role in explaining the 
observed dispersion for a given spectral type/luminosity class. 

We checked that binarity is not a problem for our spec- 
troscopic analyses. De Becker et al. (2009) concluded that 
HD 15570, HD 16691, and HD 14947 are very likely single be- 
cause the authors did not find signs of radial velocity variations 
that they could unambiguously relate to companions. HD 668 1 1 
has been extensively observed but no companions have been de- 
tected for this star, even in high-sensitivity speckle surveys for 
binarity (see e.g. Mason et al. 1998). The adaptive optics survey 
for faint companions led by Turner et al. (2008) also concluded 
that HD 210839, HD 192639 and HD 163758 are single stars. 
HD 190429 is a well-known binary but the spectrum (and the 
photometry) we used was for HD 190429A (i.e., for the primary 
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Table 1. Basic parameters of our targets. 



Star Other name Spetral Type U B V J H K E(B - V) dist. (kpc) M v 



HD 16691 




04 If 


8.620 


9.182 


8.702 


7.653 


7.444 


7.362 


0.800 


3.31±0.20 


-6-38! : 
_ 6 _23+8.n 

" -"-0.20 

-6.28+J! 


HD 66811 


<TPup 


04 I(n)fp 


0.890 


1.941 


2.210 


2.790 


2.955 


2.968 


0.040 


0.46±0.04 


HD 190429 A 




04 If 




7.201 


7.088 


6.189 


6.162 


6.150 


0.460 


2.45±0.20 


HD 15570 




04 If 


8.391 


8.796 


8.110 


6.477 


6.310 


6.158 


0.966 


2.34±0.11 


-6.73!; 


HD 14947 




04.5 If 


7.850 


8.452 


7.998 


7.037 


6.945 


6.861 


0.730 


3.00±0.20 


-6 43 +}1 " 
-643+8:15 

-0.12 

-6 39 +!i -" 


HD 210839 


A Cep 


06 I(n)fp 


4.620 


5.242 


5.050 


5.053 


4.618 


4.500 


0.513 


0.95±0.10 


HD 163758 




06.5 If 


6.458 


7.346 


7.318 


7.194 


7.163 


7.157 


0.300 


3.60±0.20 


HD 192639 




07.5 Iabf 


6.830 


7.455 


7.116 


6.300 


6.271 


6.217 


0.620 


2.00±0.20 


.631-8:i? 

u -0.23 



Note : Spectral types are taken from Sota et al. (2011), while the photometry is taken from the GOS catalog (Mafz-Apellaniz et al. 2004). 
Reddening, absolute magnitude and distances are calculated using the procedure presented in Sect. 4.1. 



only). Because the companion, HD 190429B, has a spectral type 
of 09.511, the luminosity difference is about 0.7 dex (Walborn 
& Howarth 2000), implying that its contribution to the total flux 
remains small. 

The log of the spectroscopic observations is presented in 
Table 2. 

2.2. FUV and UV data 

We extracted short wavelength, high-resolution echelle (SWP) 
IUE spectra from MAST. The SWP spectra cover the spec- 
tral range, AAl 150-2000 A, at a resolving power R = 10,000. 
Only one SWP spectrum is available for HD 14947, HD 15570, 
HD 16691, HD 163758 and HD 192639, respectively. For 
HD 190429A, we used the co-added spectrum presented in 
Bouret et al. (2005). For A Cep we selected the only IUE SWP 
high-resolution spectrum that is not saturated (see Bianchi & 
Garcia 2002). 

As part of the IUE "MEGA" campaign (Massa et al. 1995; 
Howarth et al. 1995), £ Puppis was observed over nineteen days, 
corresponding to slightly more than three rotation periods (5.2 
days). We chose to stick to the MEGA campaign to avoid any 
problems with "long-term" variations, which have not been stud- 
ied extensively. For the 107 spectra from this campaign, we con- 
structed three mean spectra - one for each "rotational" cycle. 
Because these three spectra were very similar, we averaged all 
of them to form the spectrum that is used in our analysis. 

We furthermore smoothed the spectra to a resolution of 
40 kms 1 to increase the signal-to-noise ratio. The spectra of 
all stars show many narrow lines of interstellar origin. This in- 
terstellar contamination dominates the whole UV spectrum of 
HD 15570; only strong lines such as C iv i/U548-1550 or N iv 
-11718 can be used safely for this star. 

All stars but HD 668 1 1 were observed with FUSE through 
the LWRS 30" X 30" aperture. The nominal spectral resolution 
is 20,000, or about 20 kms~', and the wavelength range ex- 
tends from 905 A to 1 187 A. The processed FUSE spectra were 
retrieved from MAST and reprocessed with the final version 
of the FUSE pipeline CalFUSE 3.2.3. Individual sub-exposures 
were co-added for each segment and then merged to form a sin- 
gle spectrum, using Lindler's FUSE-REGISTER program. We 
avoided contamination by the "worm" artifact (Sahnow et al. 
2000) by using only the LiF2A spectra on the long-wavelength 
side (/Ml 086- 1183 A) of the spectrum. Finally, the co-added 
merged spectra were smoothed to a resolution of 30 kms~' to 
enhance the signal-to-noise ratio. 

For HD 66811, we used the Copernicus spectra already dis- 



cussed in Pauldrach et al. (1994) that were originally obtained 
by Morton & Underhill (1977), to which we refer the reader for 
details concerning the observations and data reduction. 

2.3. Optical data 

For all but one of the northern stars, we used data obtained 
in November 2004 with the ELODIE spectrograph, which was 
mounted on the 193cm telescope at Observatoire de Haute- 
Provence (OHP). The resolving power is 42,000 and the wave- 
length coverage extends from 3895 A to 6815 A. The exposure 
times were chosen to ensure a signal-to-noise ratio of at least 100 
at 5200 A. The reduction of spectroscopic data acquired with 
ELODIE was performed using the standard reduction pipeline 
described in Baranne et al. (1996). Each order was normalized 
by a fit to the continuum, which was specified by fitting a smooth 
spline to carefully selected continuum windows. This step was 
then followed by a complete merging of successive orders to re- 
construct the full spectrum for each star. 

The spectrum of A Cep was collected with the NARVAL 
spectropolarimeter at Telescope Bernard Lyot (TBL) at Pic du 
Midi Observatory in December 2006. The resolving power is 
R=65,000 for a wavelength coverage 3700-10050 A. The set of 
four sub-exposures was processed using Libre ESpRIT (Donati 
et al. 1997), a fully automatic reduction package installed at TBL 
for optimal extraction of NARVAL spectra. The signal-to-noise 
ratios per 2.6 kms velocity bin range from 400 to 1100 over 
the 3800-7000 A range that we used here. 

The spectrum of HD 163758 was extracted from the UVES 
Paranal Observatory Project, a program of acquisition, reduc- 
tion, and public release of stellar spectra obtained with UVES 
at the VLT (Bagnulo et al. 2003). Complete details concerning 
these data can be found at the UVES POP URL l . The spectrum 
was normalized by fitting a spline to the continuum, which was 
defined by the same windows used to normalize the ELODIE 
spectra. 

The other southern star of this sample, namely £ Puppis, was 
observed with the FEROS spectrograph on the MPI 2.2m tele- 
scope at La Silla Observatory (prog. ID 074.D-0300(A)). The 
spectrum coverage is from 3800 to 8800 A, with a resolving 
power of 48,000. The spectrum was wavelength-calibrated and 
extracted using a modified FEROS pipeline that includes the 
modifications already described in Sana et al. (2006) and Sana 
(2009). 



1 http : //www . sc . eso . org/santiago/uvespop/ 
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Table 2. Star sample and observations logs. 



FUV UV Optical 



Star 


Spectral Type 


Data Set" 


Date 


Data Set 


Date 


Instrument 


Date 


HD 16691 


04 If 


F-E8050101 


2004-11-18 


SWP46545 


1992-12-21 


OHP/ELODIE 


2005-11-05 


HD 66811 


04 I(n)fp 


C - C044-001 


1973-02-22 


SWP15296 


1981-10-20 


ESO/FEROS 


2005-12-20 


HD 190429 A 


04If 


F-P1028401 


2000-07-18 


SWP b 




OHP/ELODIE 


2004-08-27 


HD 15570 


04If 


F-E0820101 


2005-11-08 


SWP 11234 


1981-02-04 


OHP/ELODIE 


2004-11-17 


HD 14947 


04.5 If 


F - E0820201 


2004-30-09 


SWP 10724 


1980-12-01 


OHP/ELODIE 


2004-11-20 


HD 210839 


06 I(n)fp 


F- PI 163101 


2000-07-22 


SWP52623 


1994-10-23 


TBL/NARVAL 


2006-12-13 


HD 163758 


06.5 If 


F-P1015901 


2000-08-16 


SWP02892 


1978-10-09 


ESO/UVES 


2001-07-27 


objectobjectHD 192639 


07.5 Iabf 


F-C17101010 


2002-09-04 


SWP09493 


1980-07-11 


OHP/ELODIE 


2001-08-11 



a Observatory used to obtain spectra: F FUSE; C Copernicus 
b Co-addition of several spectra; see Bouret et al. (2005) 



Some broad spectral lines like the He n /14686 line in the 
FEROS spectrum of £ Puppis or the Ha line in either ELODIE or 
FEROS spectra extend over two adjacent orders. In these cases, 
we found that adopting order shapes from adjacent orders gave 
the best continuum reduction. 

2.4. Some considerations about spectroscopic variability 

The observations used for this study have been obtained at very 
different epochs for each star (see Table 2). Therefore, we had 
to address the problem of the expected variability of the stellar 
winds, and particular how this would affect the physical quanti- 
ties that we measured through spectral fitting. 

Resonance line variations are well documented for several 
O-type stars in the FUV-UV range, including two of our tar- 
gets, namely ( Puppis and A Cep. For the other stars only one 
FUV-UV spectrum has been obtained, and hence no information 
is available on the UV variability of their wind or photosphere. 
The most prominent signature of variable wind structures in hot 
stars is the presence of blueward-migrating "discrete absorption 
components" (DACs), which are clearly visible in time-series 
IUE spectra (e.g. Howarth et al. 1995; Massa et al. 1995). These 
data have established that the variations in the wind lines are not 
chaotic fluctuations, but are instead very systematic and prob- 
ably associated with large-scale perturbations that arise from 
corotating interacting regions (Cranmer & Owocki 1996). 

In the optical range short-timescale (of about hours) stochas- 
tic spectral line variations of Ha and He n /14686 are well docu- 
mented for O-type stars (e.g. Prinja & Fullerton 1994; Eversberg 
et al. 1998; Markova et al. 2005; De Becker et al. 2009). This 
variability is commonly attributed to the small-scale structures, 
or clumps, which are predicted by time-dependent models of the 
line-driven instability (Owocki 1994; Dessart & Owocki 2003, 
2005), though pulsations might also be involved at some level 
(e.g. Markova 2002; Prinja et al. 2006). The observed line profile 
variability indicates that the winds are not smooth, but perturbed. 
Observations furthermore showed that significant variations in 
Ha start from the wind base up to velocities of a few hundreds 
of km s _1 , which corresponds to the maximum extension of the 
Ha line formation region in dense enough winds such as those 
of O supergiants. Evidently, these clumps can coexist with the 
large-scale perturbations (DACs) observed in the UV domain. 

For our modeling, we assumed that the clumping properties 
of the wind are determined by small-scale structures that origi- 
nate from the intrinsic instability of the line-driving mechanism. 
Hence, the derived clumping properties should be (almost) inde- 
pendent of time, as long as the major wind characteristics remain 



largely constant. This is supported by observations that on longer 
time scales, the wind properties of most O stars appear to be con- 
stant. Additional support comes from a study showing that for 
early-type supergiants, variations observed in Ha indicate mass- 
loss rate changes of about ± 4% with respect to the mean value 
of M (Markova et al. 2005). These changes are smaller than the 
uncertainties of our M estimates, most of which are attributable 
to uncertainties in the stellar luminosities. 

3. Model atmospheres 

We performed the analyses using model atmospheres calculated 
with the code CMFGEN (Hillier & Miller 1998). This code com- 
putes line-blanketed, non-LTE (NLTE) models by solving the 
radiative transfer and statistical equilibrium equations in the co- 
moving frame of the fluid, for a spherically symmetric outflow. 
To help facilitate the inclusion of line blanketing, super-levels 
were used. This approach allows the inclusion of many energy 
levels from ions of many different species. After convergence of 
the model, a formal solution of the radiative transfer equation 
is computed in the observer's frame (Busche & Hillier 2005), 
thus providing the synthetic spectrum for comparison to obser- 
vations. 

For the photospheric density structure we initially used the 
hydrostatic density structure for a model computed with Tlusty 
(Hubeny & Lanz 1995; Lanz & Hubeny 2003), while for the 
wind we used a standard /^-velocity law that is connected to the 
hydrostatic density structure just above the sonic point (at ap- 
proximately 15 km s~'). The mass-loss rate, density, and velocity 
are related via the continuity equation. In more advanced mod- 
els we iterated on the density structure so that the hydrostatic 
equation, 

dP/dr = -g gmv +g md (1) 
where 

Sgmv = GM/r 2 (2) 
and the radiative acceleration g ra d 

4n r 

grad = —P I XvH v dv (3) 

was satisfied. We solved for the density structure by integrat- 
ing the hydrostatic equation using the Runge-Kutta method. To 
facilitate the integration, the Rosseland LTE opacity was com- 
puted as a function of temperature and electron number density 
for the model abundances prior to the CMFGEN run. In a typ- 
ical CMFGEN model, a hydrostatic iteration is performed after 
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the first model iteration and then every n (with n ~ 8) iterations. 
Generally 3 to 4 hydrostatic iterations are needed to obtain con- 
vergence to better than 5% everywhere. 

A depth-independent microturbulent velocity is included in 
the computation of the atmospheric structure (i.e., temperature 
and level populations). A value of 15 kms~' is adopted as the 
default in our computations. 

We have computed clumped wind models with CMFGEN. 
Clumping is implemented through a parameteric, volume filling- 
factor approach, which assumes that the clumps are optically 
thin to radiation and the inter-clump medium is void. Under 
these assumptions, the clumped wind density, p(r) is related to 
the homogeneous (unclumped) wind density and the volume fill- 
ing by p(r) = pi f. The filling factor decreases exponentially 
with increasing radius (or, equivalently, with increasing veloc- 
ity): / = foo + (1 - foo) exp(-v/v c i) where v c i is the velocity at 
which clumping starts. We tuned v c i to improve the fit to some 
observed lines (Bouret et al. 2005). A limitation of this parame- 
terization is that /is monotonic whereas theoretical simulations 
by Runacres & Owocki (2002) and observational studies (Puis 
et al. 2006; Najarro et al. 2011) suggest it is non-mono tonic. 

We have accounted for shock-generated X-ray emission in 
our models 2 . An important consequence of the X-ray and EUV 
shock radiation is enhanced photoionization which results in 
"wind super-ionization" — that is the presence of highly ion- 
ized ions (such as O vi) that are not expected to be produced 
by the photospheric radiation field. These super-ions are primar- 
ily produced by the Auger process (Cassinelli & Olson 1979). 
In this process a X-ray photon causes the ejection of an in- 
ner shell electron from an ion. The ion, which is in a highly 
excited state, usually autoionizes. For CNO elements, the net 
change in the charge of the X-ray-absorbing ion is generally 
+2. Auger processes, and direct ionization by the shock radia- 
tion field, are accounted for in calculating the wind ionization. 
Only two stars in our sample have been observed with X-ray 
satellites, namely £ Puppis and A Cep with ROSAT; the X-ray 
luminosities are logL x = 32.43 and logL x = 31.92, respec- 
tively (see Oskinova 2005). The luminosity ratios, logLx/^bob 
for these stars were computed using the bolometric luminosities 
as determined from this study. For HD 190429 A, HD 14947, 
HD 15570 and HD 16691, we have adopted the X-ray luminos- 
ity of £ Puppis , while the X-ray luminosity of A Cep was chosen 
as a proxy for HD 163758 and HD 192639. 

4. Diagnosing stellar and wind parameters 

4. 1 . Stellar luminosity 

The stellar luminosity L is one of the key input parameters to 
CMFGEN. It is usually derived from the absolute magnitude, 
which requires accurate knowledge of the distance of a star. 

Because HD 15570 is a member of the cluster IC1805, we 
adopted the distance modulus of the cluster (Massey et al. 1995), 
DM = 1 1.85 ± 0.1 (which corresponds to a distance of 2.3 kpc), 
to derive My, hence L. A significant absorption is also present 
toward this cluster, which translates into a fairly high reddening 
of the stellar spectra. Remarkably, the color excess we derive 
from fitting the observed SED differs by less than 1 % from the 
E(B - V) computed from observed (B-V) and theoretical (B-V)o 
from Martins & Plez (2006). 



2 In these stars the clumps are probably primarily responsible for pro- 
ducing the observed N v resonance doublet, but in cooler stars the inter- 
clump medium is also likely to be very important (Zsargo et al. 2008) 



For all other stars, the distance is poorly known, with some- 
times fairly large differences in estimates of the distance mod- 
ulus. This leads to significant uncertainties in the luminosity, 
which in turn translates into differences in some parameters such 
as the mass and the mass-loss rate, which depend on the stellar 
luminosity through the implied value of the stellar radius (see, 
e.g., Puis et al. 1996). We adopted the following procedure: 

1. We adopted an initial absolute magnitude (and correspond- 
ing distance) for each object from the literature. We used 
My values listed in Markova et al. (2004) for HD 14947, 
HD 16691, ( Puppis, HD 192639 and A Cep, while My 
was taken from Repolust et al. (2005) for HD 190429A. For 
HD 163758, we found no My listed in the literature and we 
adopted My from the calibration by Martins & Plez (2006). 
From My, we derived L assuming a bolometric correction 
from Martins et al. (2005). We then searched for the model 
that best fits the normalized UV and optical spectra. This 
model yields the photospheric and wind parameters. We de- 
rived M by requiring that the model fits Ha. Any value of M 
and R* such that M /Rj 5 is constant provides a similar fit 
to Ha (and He n /14686). However, UV lines do not follow 
the same trend, and this allowed us to break the degeneracy 
between M and R* to some extent.. 

2. From the best-fit model, we fitted the flux calibrated IUE 
spectra and UBVJHK fluxes (for the photometry see Table 
1) to constrain the distance and the extinction. The Galactic 
reddening law of Cardelli et al. (1989) was used. We var- 
ied the distance (global scaling) and color excess E(B - 
V) (wavelength-dependent leverage) to reproduce the UV- 
optical-near-infrared spectral energy distribution (SED). 

3. The absolute magnitude My was recomputed with the dis- 
tance derived in step 2 from the fit to the observed fluxes. 
If the new value differed from the value initially assumed in 
step 1 by 0.05 dex or more, it was used to compute a new 
luminosity and a new model atmosphere. When the revised 
model atmosphere is calculated, T e g, logg, Voo, j3 and the 
abundances were kept fixed to their respective values derived 
in step 1, but M and / (the clumping filling factor) were var- 
ied to determine a new best fit to the normalized spectra. 

4. This set of three steps was repeated until we converged on a 
solution satisfying both the fit to the individual spectral wind 
features from the normalized spectra and the absolute UV + 
optical/near-IR (NIR) fluxes. 

Final values for the distances and absolute magnitudes are 
listed in Table 1 . The SEDs for each star resulting from the pro- 
cess above are presented in Fig. 1. 

4.2. Effective temperature, log g and helium abundance 

Our primary T t g diagnostics in the optical were the ratio of He i 
/J4471 to He n /14542 equivalent widths, as well as line profile 
fitting. We used the other lines of He i and He n to refine the 
r e ff determination. The optical helium lines were also used to 
constrain the relative He to H abundance. In addition to these 
classical helium diagnostics, we also used metal lines for esti- 
mating the effective temperature. This is particularly relevant for 
the earliest stars of the sample, which have weak He i lines, and 
for which the validity of using the helium ionization balance to 
determine T e g accurately is questionable. In the UV the Fe iv 
to Fe v and Fe v to Fe vi ionization ratios provide a sensitive 
temperature diagnostic. We found consistent estimates with the 
helium ion ratios within one thousand Kelvin for most objects, 
a slightly higher difference for £ Puppis and A Cep is related to 
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Fig. 1. Synthetic SEDs (in red) compared to flux calibrated + UBVJHK photometry for the target stars (in black). The distance, E(B-V) and 
luminosity were iterated to reach agreement between models and observations (see Sect. 4.1). For plotting purpose, the fluxes were scaled by a 
factor 10 14 for all stars but HD 6681 1 and HD 210839, where factors 10 9 and 10 13 were used, respectively. 



the higher v sin i of these objects. When possible ( HD 190429 A 
and HD 16691), we also used the C iv A 1169 to C m ^1176 line 
ratio which has also been shown to provide a useful temperature 
diagnostic by Heap et al. (2006). 

In some cases, effective temperature determinations based on 
the ionization ratio of optical helium lines are inconsistent with 
some temperature-sensitive metal lines, such as N v AAA6Q5- 
4620. The presence of these lines in the observed spectra points 
to an enhancement in T e g of approximately 2000 K compared to 
the values derived from the helium lines. The relative strength 
of the N v doublet, which is strongly dependent on the temper- 
ature, has a relatively weak dependence on gravity, while the 
absolute strength was better reproduced by tuning the nitrogen 
abundance. Nevertheless, we could not achieve a good fit to these 
lines without degrading the agreement with the rest of the opti- 
cal and UV spectrum. The values of T e g and nitrogen abundance 
needed to fit these lines should be regarded as upper limits. 

The Stark-broadened wings of Hy provided the primary con- 
straint on \ogg (other Balmer lines were used as secondary in- 
dicators). The typical uncertainty in logg is 0.1 dex. Despite the 
high S/N of our spectra it is difficult to do better than this — 
logg has to be determined simultaneously with T e ff, and there 
are slight inconsistencies with different lines (because of rectifi- 
cation errors and residual uncertainties with the atomic physics). 

4.3. Micro-turbulence 

The formal solution of the equation of transfer was used with a 
radially dependent turbulent velocity to calculate the emergent 
spectrum (see also Hillier & Miller 1998). In this case, the mi- 



croturbulent velocity follows the relation £ t (r) = (£™ flx - 
where and £" lax are the minimum and maximum mi- 
croturbulent velocities. The microturbulence at the photosphere 
is chosen to best reproduce the strength of the S v A 1502 and iron 
lines (Fe iv and Fe v) in the IUE spectra. The adopted v tul b has 
little effect on the strength of most H i lines. On the other hand, 
we found some significant changes for important T e g diagnos- 
tics such as He i AAAll and some He n lines (see also Villamariz 
& Herrero 2000). We confirm findings by Hillier et al. (2003) 
that He n /15411 is also sensitive to the microturbulent velocity, 
while He n A<\5A\ or He n /14686 are not. 

In the outer region of the wind, the turbulence is adopted to 
best fit the shape and slope of the blue side of the absorption 
component of the C iv AA 1548- 1551 P Cygni profile. 

4.4. Surface abundance 

Accurate CNO abundances provide key constraints on both the 
evolutionary state and the influence of rotation, and should 
preferably be determined from photospheric lines of these ele- 
ments. 

We adopted the solar abundances from Asplund et al. (2005) 
as a reference. The new values better agree with surface abun- 
dances of B-type stars in the solar neighborhood (see Cunha 
et al. 2006; Przybilla et al. 2008; Simon-Diaz 2010; Nieva & 
Simon-Diaz 2011). Some caution may be required when com- 
paring these baseline abundances to those we derives since only 
one star ( £ Puppis) is located within 1 kpc, which is the limit 
usually set for solar neighborhood abundance studies (but see 
the recent work by Nieva & Przybilla 2012). 
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4.4.1. Carbon abundance 

Most carbon lines in the FUV/UV range suffer from wind con- 
tamination and we chose not to use them as prime indicators for 
carbon abundances. For stars with T e g < 35,000 K, C in ,11247 
appears as an emission line with a clear response to abundance 
variations. This is the only diagnostic for carbon abundance in 
the UV for these stars. 

In the optical, available lines depend on the temperature 
range (except for C iv ,1,15801-5812 which are present for all 
stars). In a few stars, C in /L14647-50-52 lines are clearly de- 
tected and were used. For stars cooler than 38,000 K, C m ,15696 
was used. For stars with T eS < 35,000 K, C in ,14068 is present 
as a non-saturated absorption line and was also used to measure 
the carbon abundance. Finally, we checked for consistency with 
the other C in diagnostics listed above using the C in emission 
lines at ,16721 , C in /U6727-6730. 

4.4.2. Nitrogen abundance 

In the UV, the photospheric lines of N in /M1748-1752 were 
used. However, these features saturate when the nitrogen abun- 
dance exceeds about five times the solar value, at which point 
they are no longer useful abundance diagnostics. 

In the optical, N hi ,14379, N ra ,1,14634-4640-42, N m 
/L14510-4515, N in ,15320-5327, N iv,U5204-5205, N iv,16380, 
N v ,14944 and N iv ,16212-6220 are available. We furthemore 
used the emission lines of N v ,1,16716-6718 to check for consis- 
tency only, because these lines are often blended with C m lines 
(/16721, ,1,16727-6730). 

We found that the N iv ,15205 line provided the most reliable 
abundances in optical spectra of 04 - 05 supergiants. It is sensi- 
tive to abundance changes, with no sign of saturation within the 
limits we tried. While there may be problems with the absolute 
calibration of abundances based on this line, relative abundances 
should be much more accurate. 

4.4.3. Oxygen abundance 

The O m ,15592 line provides the best constraint on the oxygen 
abundance in optical spectra. Indeed, this line is not contami- 
nated by wind contribution and is observed as an absorption line 
(though often weak), in a region of the spectrum where the con- 
tinuum is well defined. We also used the O m triplet (AA1 150- 
1 154 A) in the FUSE spectrum, when fast rotation does not blend 
them with other metal lines. UV lines (O iv ,1,11338-1343 and 
O v ,11371) indicate a (sometime significant) contribution from 
the wind and were not used. 

4.4.4. Formation of some specific lines 

Before proceeding further it is worth making some general com- 
ments about photospheric lines, which sometimes appear in 
emission in O stars. If we assume the Eddington-Barbier relation 
holds, these transitions are in emission because the source func- 
tion is larger than the neighboring continuum source function. 
This can occur because the population of the upper level is en- 
hanced, or because the population of the lower level is depleted. 
Many of the transitions occur between the n=3 levels of CNO 
ions (n=4 for Si iv), and involve one, or more levels, coupled 
to low-lying levels by transitions in the EUV (i.e., shortwards 
of 900 A). Depending on the EUV radiation field, these transi- 
tions can drain or populate a particular level. These transitions 



can also be influenced by overlapping lines from other species 
(e.g. Najarro et al. 2006), and by the adopted abundances (which 
directly alters the optical depth of the line). In CNO elements, 
dielectronic recombination may play a role by selectively en- 
hancing the populations of some levels relative to others (Rivero 
Gonzalez etal. 2011). 

In WN stars, Hillier (1988) found that the strong N iv emis- 
sion lines (4058 A, 3478-3485 A, 7013-7129 A) are influenced 
by continuum fluorescence and dielectronic recombination. The 
C iv ,1,15801-5812 doublet is also strongly influenced by con- 
tinuum fluorescence, which is pumped by a transition at 312 A. 
On the other hand, in WC stars, with their large carbon mass- 
fraction (> 0.1), radiative and dielectronic recombination domi- 
nate the formation of optical carbon lines (Hillier 1989). For O 
stars the situation is likely to be more complicated (the densi- 
ties are higher and more processes are likely to be important) 
because LTE must be recovered close to the line formation re- 
gion. Fits to these lines significantly improved when we com- 
puted models with a much larger number of atomic elements, 
which were originally included to improve the calculation of the 
line radiative force. We found that the change induced in the to- 
tal EUV line-blanketing yielded C iv ,1,15801-5812 profiles that 
better agreed with observations. 

C in ,1,14647-50-52 lines form very near the stellar surface 
and are marginally sensitive to wind parameters. In addition, 
the lower level of C m ,1,14647-50-52 (ls 2 2s3s - ls 2 2s3p) cor- 
responds to the upper level of C in transitions around 538 A 
(ls 2 2s2p - ls 2 2s3s). Around this wavelength, Fe iv lines are 
present and overlap with the C in line. These Fe iv lines drain 
and depopulate the lower level of the C in AA4647 -50 transition 
(see also Martins & Hillier 2012, in preparation). We had to re- 
duce significantly the oscillator strengths of Fe iv lines around 
538 A to obtain satisfactory fits to these C m lines. This sensi- 
tivity to atomic model assumptions must be kept in mind when 
using these lines for abundance determinations. 

C in ,15696 requires many energy levels in our model atom 
to form in emission as observed. However, in some cases we 
could not reach a satisfactory fit of the lines, indicating that some 
mechanism or process that is affecting the formation of this line 
has not been correctly accounted for in our modeling. To our 
knowledge, systematic studies of the formation of C m emission 
lines have been performed for WC stars only. We note first that 
C in ,15696 behaves differently from other C m emission lines, 
even though it varies smoothly along the WC spectral sequence, 
(e.g. Torres et al. 1986). C m lines such as ,12296 and C in ,16741 
are strong in WC4 stars, while ,15696 is very weak or absent. On 
the other hand, in late WC stars, ,15696 is the strongest optical 
C ni transition. This behavior occurs because the upper level of 
,15696 (2s 3d l D) can also decay to the 2s 2p 'P° via a transi- 
tion at 574 A — only when the transition is optically thick will 
,15696 be strongly in emission (A(,1574)/A(,15696)=150, Storey 
[private communication]) (Hillier 1989). The lower level is cou- 
pled to the 2p 2 'D state by a transition at 884A, highlighting the 
importance of EUV transitions on the source function for the 
,15696 transition. Note that ,14058 is the N iv isoelectronic equiv- 
alent transition to C in ,15696. 

Mihalas et al. (1972) showed that the N m ,1,14643-4640- 
4642 feature can be driven into emission by dielectronic recom- 
bination, and highlighted the importance of two-electron transi- 
tion for draining the 3p state. They also noted that the Swings 
mechanism (i.e., continuum fluorescence) could be important in 
stars with winds. Very recently, Rivero Gonzalez et al. (2011) 
showed that in stars with stellar winds, the Swings mechanism 
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is indeed more important, although they also note that in O stars 
cooler than 35,000 K the interaction between the O in and N m 
resonance lines could also be important. 

N iv ,1/16212-6220 (2s4 3 S -2s4p 3 P°) appears as an emis- 
sion feature in the earlier supergiants. Our original models failed 
to reproduce this line, which was traced to a problem with the 
theoretical model atom. In the atomic data calculations the lev- 
els designations of the 2s 4s 3 S° with the 2p 3p 3 P° level (as pre- 
viously noted by Laughlin 1990; Allard et al. 1991) are inverted 
relative to the spectroscopic designations. When we rectified 
this inversion, emission in N iv A62 12-6220 was correctly re- 
produced. 

In our models for some of the 04 supergiants, the N iv /14058 
line shows up as an intense emission feature, which is not ob- 
served. The N v 4605-4620 lines are present as photospheric 
absorptions in 04 supergiant spectra only. Since they are also 
very sensitive to effective temperature we did not use them to 
constrain the nitrogen abundance. Walborn (2001) describes an 
apparent connection between the N v ,M4605-4620 and the N iv 
/14058 lines: the former appear in absorption in early-type super- 
giants, while the latter are seen in emission of "similar inverse 
intensity". Although the presence and relative intensity of these 
lines is clearly defined in Of supergiants, there is no obvious 
physical connection that would explain their observed propor- 
tionality (other than they both depend on T e «). 

4.5. Wind parameters 

The wind terminal velocities, Voo, were estimated from the blue- 
ward extension of the absorption part of UV P-Cygni profiles, 
which occurs up to v^ + v max , where v max is the maximum mi- 
croturbulent velocity described above. Fits of the UV P-Cygni 
profiles using the above relation for microturbulent velocity al- 
lows a direct determination of Voo. The typical uncertainty in our 
determination of v^ is 100 km s _1 (depending on the maximum 
microturbulent velocity we adopt). 

Mass-loss rates are derived from the simultaneous analy- 
sis of strong UV lines such as Civ/11169 + Cni/11176, Nv 
A1240, C iv /U1548, 1551, Si iv A/11394, 1403, O v A1371 and 
N iV/11718. In the optical the strength of Ha and He n /14686 is 
sensitive to M, while the shape of their emission profile is sensi- 
tive to the j3 exponent of the wind velocity law. 

Following Bouret et al. (2005), dumping-related quantities 
(M, foo and v d ) were derived from P v M\\\%, 1 128, O v A1371, 
and N iV/11718 in the FUV/UV domain, to which we added 
clumping sensitive lines from the optical, primarily Ha and He n 
/14686. We emphasize here that some photospheric lines in the 
UV and optical are also sensitive to the adopted filling factor 
(and scaled M). For photospheric H i and He lines for instance, 
this is essentially caused by a weaker wind contribution (emis- 
sion) in clumped models, which produces deeper absorption 
compared to smooth-wind models. For a more detailed discus- 
sion concerning the sensitivity of line to the wind and clumping, 
see Hillier et al. (2003). 



5. Results of the spectroscopic analysis 

The synthetic spectra were convolved with a representative in- 
strumental profile and a rotational profile calculated for a trial 
value of the projected rotational velocity (v sin i) to permit di- 
rect comparison with observations. We varied vsin i until we 
achieved a good fit to the observed photospheric profiles. 



As an independent check, we used the Fourier transform 
method described in Simon-Diaz & Herrero (2007, and refer- 
ence therein) to determine v sin from the observed spectra. 

To apply the Fourier transform method, we selected lines 
with as weak a wind contribution and as symmetric a profile 
as possible. We identified several isolated photospheric lines 
(which included N v X14604 - 4620 for HD 190429 A and 
f Puppis, O in /15592 for cooler stars), which we used to de- 
rive individual values that we averaged to yield the values listed 
in Table 3. For the fast rotators, namely £ Puppis and A Cep, 
v sin i from the two method agree better than within 10 km s _1 , 
which is well within the typical error bars of each method. For 
the other moderately rotating stars, the difference is smaller than 
20 km s _1 , which is slightly smaller than the difference found for 
other O-type supergiants by Simon-Diaz & Herrero (2007). The 
final values used in our models are listed in Table 3. 

The quality of the fits was then additionally improved by 
convolving the rotationally broadened profiles with a Gaussian 
of FWHM 

Vmac to account for (isotropic) macro -turbulence. 
More complex description of the macro-turbulence exists, e.g. 
the so-called radial-tangential broadening. In the specific case of 
massive stars, Aerts et al. (2009) recently interpreted the miss- 
ing line broadening in terms of the collective effect of numerous, 
low-amplitude, non-radial gravity-mode oscillations. These au- 
thors emphasized that using the Gaussian formulation of macro- 
turbulence to fit line profiles could lead to significant underesti- 
mation of the actual rotational velocity and warned that appro- 
priate expression for the pulsational velocities should be used 
instead. Therefore, the values quoted in Table 3 should be re- 
garded as lower limits, but we do not expect significant modifi- 
cations of our conclusions concerning surface abundances. First, 
rotation and macro-turbulence in the adopted description are not 
expected to affect line equivalent widths, only the line profiles 
(contrary to micro-turbulence). Second, differences between ro- 
tational broadening from line fitting, which should be regarded 
as an upper limit, or from Fourier transform are quite small (less 
than 20 kms -1 , see above). Third, in the total broadening we 
considered, the rotational term is always dominant with respect 
to macro-turbulence. 

Given the complexity of the fitting process over the FUV 
to optical spectral range, we chose not to derive formal statis- 
tical uncertainties, nor did we try to estimate the correlation of 
errors between parameters. Instead, we varied the different pa- 
rameters until we obtained the solution that provided the best 
fit by eye to the observations, with more weight being given to 
primary diagnostics (cf Sect. 4). Then, for a more quantitative 
estimate, we varied each parameter independently around this 
solution and computed the residual (i.e. the difference) between 
the observed and the synthetic spectrum. We imposed the con- 
dition that this residual must remain within ±10% and adopted 
the range of values that fulfilled this condition as error bars on 
each parameter. The stellar and wind parameters we derived are 
gathered in Table 3. We quote all chemical abundances by num- 
ber and mass fraction, except for He i, which is given by number 
relative to hydrogen (see Table 4). As for the uncertainty in the 
CNO abundance measurements, error bars are based on the fit 
quality of the lines listed above, as estimated from a classical^ 2 
procedure. Our conclusions for individual objects are summa- 
rized in the remainder of the section. The best-model fits to the 
data are presented in Appendix A. 
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5.1. HD 1 90429 A -04 If 

HD 190429A was studied extensively in the UV in Bouret 
et al. (2005). We recomputed their best-fit model, using the new 
Asplund et al. (2005) abundances. We scaled the CNO as well 
as P and S abundances to the same relative values as in Bouret 
et al. (2005). The best fit to the P v AAl 1 18-1 128 resonance dou- 
blet is obtained for P/P Q = 0.7, further confirming the apparent 
depletion in phosphorus of this star. Only a limited tuning of the 
model parameters was needed to improve the fit to the optical 
lines, which demonstrates the reliability of our FUV/UV anal- 
ysis. The best fit to FUV-UV and optical spectra is presented 
in Appendix A. Note however that the presence of strong lines 
of N v /1/14605-4620 would require significantly higher tempera- 
ture to be reproduced (r e ff 43,000 K), which would then make 
the fit worse everywhere else in the spectrum (see also comments 
on Sect. 4.2). 

The C m /O4647-50-52 lines are weak, especially compared 
to the strong neighboring N m ,1,14634-4640-4642 lines. This 
was recognized (e.g. by Walborn et al. 2010) as due to an 
abundance effect, as CNO processed material is expected at 
the surface in supergiants. We found that the atmosphere of 
HD 190429 A is highly depleted in carbon (X(C) = 0.009% by 
mass), and is enriched in nitrogen (X(N) = 0.73%). Although de- 
rived from fewer lines, the oxygen abundance corresponds to a 
mass-fraction X(O) = 0.067%. Together with the He enrichment 
y = 0.2, this is qualitatively consistent with the general trend 
predicted by evolutionary models with rotation. 

5.2. £ Puppis - 04 1(n)fp 

Although a model with T e tf = 39,000 K provides an excellent fit 
to helium lines, raising T t g up to 41,000 K better agrees with 
the observed N v ,04605-4620 lines, which are otherwise too 
weak. Deriving the actual value proved to be elusive, but we ex- 
pect that T eB = 40,000 ±1000 K is a realistic estimate of the 
effective temperature of £ Puppis. The wind properties agree 
very well with those presented in Najarro et al. (2011) based 
on the modeling of the near-infrared spectrum of £ Puppis. Only 
clumped-wind models matched the whole set of wind-sensitive 
lines. The best fit to P v/Lll 1 18-1 128 lines was nevertheless ob- 
tained only when using sub-solar phosphorus abundance, P/P Q 
= 0.5. We needed to start clumping at v c i = 100 kms~', signifi- 
cantly above the sonic point to obtain a reasonable fit of the cen- 
tral trough of Ha. This indicates that more absorption is needed 
in the line formation region, which is missing if clumping starts 
right above the photosphere. It is remarkable that this behav- 
ior is observed on the other fast rotator of our sample A Cep, 
as well as in HD 16691, whose emission lines present strong 
similarities with those of £ Puppis (see below): N m /04634- 
4640-4642 presents a box-shaped profile, while He n ,14686 is 
clearly double-peaked. Based on preliminary models using a re- 
alistic 2D treatment of rotation as described in Busche & Hillier 
(2005) and Zsargo et al. (2006), we concluded that the shapes 
of these emission lines are most easily interpreted in terms of 
an effect of high rotational velocity. These conclusions will be 
presented in a forthcoming paper. 

Despite the high v sin i of this star, many lines of CNO ions 
are present in the spectra to derive reliable estimates of the sur- 
face abundances (see Fig. A.4). One such line worth mentioning 
is the N iv ,07103-7129 multiplet, which is the triplet equiv- 
alent of N iv ,14058. Our FEROS spectrum of £ Puppis shows 
an emission feature at these wavelengths whose strength is very 
sensitive to abundance variations and is well-matched in our fi- 



nal model. The total of CNO abundances is clearly super-solar 
in £ Puppis (see Table 4). It has been argued that this could be 
the result of a binary/merger evolution (Pauldrach et al. 2012), 
in line with the runaway status often invoked for this star (Sahu 
& Blaauw 1993). We stress, however, that in Pauldrach et al. 
(2012), this conclusion is based on abundances for several ele- 
ments, which are very different from what we derive. Their high 
C abundance for instance seems impossible to reconcile with 
the observed optical spectrum when modeled with CMFGEN. 
Detailed comparisons between the two analyses are needed to 
elucidate this problem. 

We also note that the luminosity we adopted differs signif- 
icantly from those quoted in Najarro et al. (2011) or Pauldrach 
et al. (2012). These two studies chose the higher value, often re- 
lated to its runaway status, while we adopted a lower value for 
the distance, compatible with the Hipparcos distance quoted for 
f Puppis. On the other hand, Schilbach & Roser (2008) found 
f Puppis to be even closer (300 pc) than given by Hipparcos 
(430 pc). These authors suggested that the star originally was in 
the open cluster Trumpler 10, which it left about 2.5 Myrs ago. 
A direct consequence of this lower distance would be to revise 
the luminosity downward to log L/L Q = 5.42. The stellar mass 
would go down accordingly (« 17 M ), which starkly contra- 
dicts the usual parameters of an 04 supergiant. 

5.3. HD 16691 -04If 

A model with T eff = 41,000 K ± 1000 K best fits the He lines 
as well as N v ,04605-4620 in HD 16691. Carbon lines in the 
optical are weak, which suggests a high depletion of this ele- 
ment (but see Sect. 4.4 for comments about the sensitivity of C iv 
,05801-5812 doublet to the treatment of UV blanketing). This 
is also supported by the C iv Al 169 and C m ,11 176 lines, which 
show up as weak absorptions. Nitrogen lines are strong, indicat- 
ing that this element is altogether highly enriched. We note that 
despite several tests we performed, N iv ,14058 is predicted to be 
strongly in emission compared to the observed profile. Fitting 
this line would require a much lower nitrogen abundance, incon- 
sistent with the rest of the nitrogen lines. In accordance with its 
Of/WN transition status (Conti & Bohannan 1989), HD 16691 
furthermore presents a clear enrichment in He. 

Only clumped models agree well with the observed wind 
sensitive lines. The best fit to P vAAl 118-1 128 lines was ob- 
tained using sub-solar phosphorus abundance P/P Q = 0.5. 

UV and optical spectra of HD 16691 and HD 190429A show 
strong similarities, reminiscent of those seen in their IR TT-band 
spectra (e.g. the subordinate He n 2.189yum in emission). Conti 
et al. (1995) noted that both stars would have been classified 
as WN stars had they been observed only in the TT-band. These 
spectral similarities strongly suggest a common evolutionary sta- 
tus for these objects. 

HD 16691 and £ Puppis also share some striking similarities 
(see Table 3), especially regarding the emission lines of their op- 
tical spectra, such as N m ,04634-4640-4642 and He n ,14686. 
We argue that a high rotation rate is able to explain the observed 
emission profiles of HD 16691, and there is no need to invoke 
an extended corotating region that would be caused by a pu- 
tative magnetic field, as proposed by De Becker et al. (2009). 
Furthermore, we needed to start clumping at v c i = 130 kms~' 
in our models to best reproduce the observed spectra. This be- 
havior is again very reminiscent of that of the other fast rotators 
of our sample (see below) and might indicate that the wind of 
HD 16691 is indeed rotating faster than suggested from the pho- 
tospheric lines. 
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5.4. HD 15570-04 If 

HD 15570 is one of the brightest early-type stars in the young 
open cluster IC 1805. The FUV-UV spectra are heavily masked 
by interstellar absorption lines (mostly H2 in the FUSE spec- 
trum), consistent with the large E(B - V) we derive (see Sect 
4.1). The microturbulent velocity, which is usually best deter- 
mined from S v A 1502 and the iron line strengths (Fe iv and 
Fe v) could not be accurately measured and we used optical lines 
only, to derive a value g t = 15 km s . 

We derived r eff = 38,000 K, log L/L Q = 5.94 for HD 15570, 
which implies that it is somewhat cooler than expected for 
its spectral type, although at the expected luminosity accord- 
ing to Martins et al. (2005). Usually the observed differ- 
ences/dispersion in this calibration can account for a large part of 
the uncertainty in the luminosity, i.e., the distance to this object. 
This should not be the case for HD 15570 because the distance 
modulus IC1805 cluster is well determined. 

More puzzling, the He lines of HD 15570 are weak for its 
spectral type. Since the photospheric and wind parameters of 
HD 15570 (see Table 3) are very close to those of HD 190429A, 
we compared the optical spectra of both stars. The unusually 
weak helium lines of HD 15570 are clearly apparent in Fig. 2, 
where He lines are stronger and broader in HD 190429 A (v sin i 
for HD 190429A is higher than for HD 15570). 



1.030 




0.92 ■ 

0.90 I , , , , I I , , 

5860 5865 5870 5875 5880 5885 6670 6680 6690 6700 
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Fig. 2. Helium line profiles from the optical spectra of HD 190429 A 
(black) and HD 15570 (red). HD 190429A has vsin i = 150 kms" 1 , 
while it is 97 km s _I for HD 15570. 

The high-resolution spectrum of HD 15570 presented in De 
Becker et al. (2006) 3 was kindly provided to us by G. Rauw for 
comparison to our spectrum. Both spectra are remarkably simi- 
lar throughout the common spectral range (4000 - 6800 A), in- 
cluding the weak and narrow(er) helium lines with respect to the 
star's photospheric properties. This confirms that our finding is 
not an observational artifact caused by incorrect data reduction 
or spectrum rectification. Spectra retrieved from the ELODIE 
archives 4 , which cover the epoch from 1996 to 2005 and have 

3 Obtained with the echelle spectrograph at the 2.1 m telescope at 
the Observatorio Astronomico Nacional of San Pedro Martir (SPM) in 
Mexico 

4 http://atlas.obs-hp.fr/elodie/ 



a signal-to-noise ratio of 100 or higher, also display the same 
narrower, weaker helium lines. This confirms that the weakness 
of He lines is not time-dependent, at least on this timescale. 

Quantitatively, we found that y(He/H) = 0.07 (by number) is 
needed to best fit the helium lines in the optical domain. Such a 
low He abundance is not supported by stellar evolution consider- 
ations, but we could not find any problems with the observations 
that would affect the analysis of He lines. It is unlikely that this 
lower helium abundance is an indication of the initial compo- 
sition of the molecular cloud from which HD 15570 formed. 
Martins et al. (2005) found y = 0.1 for HD 15629, a main- 
sequence 05V((f)) star in IC 1805, which suggests that the initial 
helium content of this cluster is normal. In the rest of this work, 
we use y = 0. 1, which should be considered as an upper limit on 
the helium abundance. 

Only models strongly depleted in carbon match the carbon 
lines, while all nitrogen lines in the optical suggest an enrich- 
ment in nitrogen. As for HD 16691, N iv/14054 appears as a 
strong emission in our models, which we could not reconcile 
with the observed profile (a very weak emission). The O in 
/15592 line shows that oxygen is also depleted. Qualitatively, 
the CNO abundances are consistent with the evolved nature of 
HD 15570, as inferred from its position in the H-R diagram. 

The best agreement for wind-sensitive lines was obtained 
with clumped models, still requiring a sub-solar phosphorus 
abundance to match P vAAl 1 18-1 128 (see Table 4). 

5.5. HD 14947 - 05 If 

The best fit was obatined for r eff = 37000 ± 1000 K, which is 
1500 K cooler than expected from the T e g - spectral type cal- 
ibration by Martins et al. (2005). This T e ff is low enough that 
the C in /15696 line is present in the spectrum and can be used 
to measure the carbon abundance, which is mildly depleted. The 
adopted value also gives a consistent fit to C in AA4641 -50-52 
and C iv /L15801-5812. We obtained a consistent fit to several 
nitrogen lines, including N m /U4634-4640-42, N m AA451Q- 
4515, N IV/L15204-5205 and N iV/16380 by adopting a surface 
enrichment comparable to those we measured in other stars of 
this sample. Similarly, O in FUV triplet (AA1 150-1 154) and O m 
/15592 indicate a depletion in oxygen. Qualitatively, the CNO 
patterns are compatible with both the evolutionary status (indi- 
cated by the position of HD 14947 in the H-R diagram), and the 
relatively minor He enhancement (y = 0.12). 

The UV and optical wind lines could only be fitted with very 
highly clumped models. The best fit was obtained for a model 
with /oo = 0.03, implying a reduction of a factor of six for the 
mass-loss rate compared to the homogeneous model. This is 
by far the highest reduction of the whole sample. This highly 
clumped model is dictated by the P v resonance doublet, for 
which we could achieve a very good fit by adopting a normal so- 
lar abundance for phosphorus. Using this low filling factor also 
maintains a good fit to all other wind/clumping-sensitive lines 
(see Sect. 4). At the temperature we derive for HD 14947, the 
effect of introducing clumping is to strongly recombine P +4 ion 
to P +3 and the former ion is no longer dominant. The opacity in 
the line is greatly reduced and thus there is no additional need to 
vary the abundance to improve the fit. 

5.6. A Cep - 06 1(n)fp 

This star is another famous rapid rotator, with conspicuous spec- 
tral variability (e.g. Howarth & Prinja 1996; Kaper et al. 1997). 
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The best fit we could achieve to the FUV/UV and optical spectra 
was for T e ff = 36,000 K, logg = 3.5. Despite the high rotation 
rate, several lines are available for accurate abundance determi- 
nations. We found that A Cep is only moderately enriched in 
nitrogen, and depleted in carbon and oxygen (Table 4), which 
is quite surprising because more extreme patterns would be ex- 
pected from rotational mixing in such a fast rotator. We also 
found that y = 0.12, in line with this moderate chemical evo- 
lution. 

The introduction of clumping in the model dramatically im- 
proved the fit to the spectrum, including Ha. Here again, we 
needed a lower-than-solar phosphorus abundance to achieve the 
fit to the Pv FUV lines presented in Fig. A. 12. As discussed 
in Sect. 6.4, this indicates that the model optical depth for the 
resonance doublet is too high for solar abundance, which likely 
implies that either the ionization balance of phosphorus is com- 
puted incorrectly by CMFGEN, or that its treatment of clumping 
is incomplete. 

The most severe disagreement between the model and the 
observations is for He n A4686 - with or without clumps we fail 
to reproduce the double-peaked profile, which is typical of a high 
rotation rate. In this context it is no surprise that the N m AA4634- 
4640-4642 lines present the same profile as already observed in 
f Puppis and HD 16691. 

5.7. HD 1 63758 - 06.5Iaf 

A low effective temperature T e g = 32,500 - 33,000 K would 
provide a very good fit to the forest of photospheric iron ions 
lines in the FUV-UV domain (see also Bianchi & Garcia 2002). 
However, such a cool effective temperature is ruled out from the 
analysis of the optical spectrum and T e g = 34,500 ± 1000 K is 
preferred to fit the He i and He n lines in the optical and the FUV- 
UV spectrum as well. Although significantly higher than the UV- 
based T e s, it is also more compatible with the spectral type of the 
star (Martins et al. 2005). We found no documented variability 
that would explain the contradictory conclusions reached from 
the FUSE - IUE spectra on one hand and optical spectra on the 
other hand. 

The carbon lines in HD 163758 are stronger than in other O 
supergiants of similar spectral type (Gomez & Niemela 1987a). 
This is the case for lines in emission (e.g. C in A1247, C m 
AA4647-50-52 and C m ,15696) or in absorption (e.g. C m A4068 
or C iv AA5801-5812). All these suggest that this star is carbon- 
rich. Our analysis shows that HD 163758 is indeed the most 
carbon-rich star of our sample, with an abundance of one and 
a half times the solar carbon abundance. While this abundance 
allows model spectra to match most photospheric C in and C iv 
lines, it still does not reproduce the strong C m A5696 emission 
line. 

HD 163758 is also one of the most nitrogen-rich stars star 
of the sample, as indicated by the strong emission lines of N m 
AA4634-4640-4642 (the N m AA1748, 1752 absorption lines are 
saturated). Because nitrogen is produced at the expense of car- 
bon, and assuming mixing processes are unchanged when super- 
solar element abundance are present, this suggests that the ini- 
tial carbon abundance of this object was very high. On the other 
hand, the oxygen abundance is sub-solar, as indicated by the very 
good fits to the photospheric O m lines in the FUV (AA1 150, 54) 
and optical (A5592) spectrum.This seems at odds with the con- 
clusion that the initial CNO mixture of this star was super-solar. 

Whatever parameters we chose to fit the photospheric fea- 
tures, the whole set of wind lines could only be matched with a 
clumped models. The improvement is especially striking for the 



Si iv AA1394-1403 resonance doublet; this suggests that for this 
range of effective temperature, Si iv in addition to P v, can be 
used to probe clumping, thus extending the conclusions reached 
by Prinja & Massa (2010) for B-type supergiants up to signifi- 
cantly higher temperatures (or different ionization fractions). 

5.8. HD 192639 -07.5Iabf 

The best fit was obtained for r e fj = 33,500 ±1000 K. In constrast 
to HD 163758, the analysis of both FUV-UV and optical spectra 
yields consistent results for HD 192639. Given the strong line 
variability of this star, especially in the He n /14686 line (e.g. 
Rauw et al. 2001), and the long time span separating the FUV, 
UV and optical spectra, this is quite remarkable. 

Our analysis indicates that HD 192639 is moderately de- 
pleted in carbon and oxygen, although its nitrogen enrichment 
is relatively stronger. Whatever carbon abundance we used, we 
could not reproduce the C m /15696 or C iv AA5801-5812 lines 
with the correct intensity (see Fig. A. 16). Note that the latter 
would be better reproduced using higher T e g, which is ruled out 
from the fit of He lines. The carbon abundance we quote is the 
upper limit above which C in A4068 and C m /L14647-50-52 be- 
come stronger than observed. The nitrogen abundance is con- 
strained from the strong N m lines in the optical at A/14510-4515 
and N m AA4634-4640-42. The N in AAl 748- 1752 lines are satu- 
rated and were not used. This abundance is consistent with other 
N in diagnostics, including N in A5320-5327 in emission, but fail 
to reproduce the weak absorption of N ivAA5 204-5205. Oxygen 
lines are altogether very well defined and responsive to varia- 
tions of abundances, either in the UV or in the optical. 

Finally, clumped models are needed to reproduce the wind 
lines, together with a factor of two depletion of the phosphorus 
abundance, compared to solar, to fit the P v FUV lines. As for 
HD 163758, the improvement in the fit resulting from clumped 
models is striking for the Si iv AA1394-1403 doublet. 

6. Discussion 

6.1. Evolutionary status 

We constructed the H-R diagram of the sample using effec- 
tive temperatures and luminosities derived from the modeling, 
as well as evolutionary tracks and isochrones from Meynet & 
Maeder (2003) (see Fig. 3), with an initial rotational velocity of 
300 kms -1 and solar metallicity. The most remarkable feature 
of this diagram is how well it distinguishes the different spectral 
types in terms of masses and ages. The 04 supergiants cluster 
tightly around the 60 M Q track, 3 Myr isochrone. Later spectral 
types probe lower masses or more evolved regions of this dia- 
gram, but in any case, the stars are more massive than 40 M 
and younger than 5 Myr. 

Alternately, the position of the 04 supergiants in this dia- 
gram might be explained by assuming that they are in fact much 
older, and located on the 60 M Q track in the range 4.3 Myr to 
4.4 Myr. However, the predicted properties of this track in this 
age range, like the hydrogen mass fraction for instance, are in- 
compatible with the observed properties of the 04 supergiants. 
Indeed, we find X(H) = 55% or higher while X(U) * 4% or 
lower should be expected for a star with initial mass 60 M Q after 
4.3 Myr. This definitively rules out this possibility. For the same 
reason, the 04 stars of this sample could not be stars with initial 
masses slightly above 85 M Q , which would be around 3.7 Myr 
old, because the expected X(H) is 5% or less. The same comment 
holds for the later-type stars though for different initial masses. 
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Table 3. Fundamental parameters and wind properties for the stars in our sample. 



Star 


r cff iogg c log^ 


log(M) Voo 


P foo v a 


v sin 


*mac 




M' 
'"evol 


M m " 


Age 


HD 


[kK] 


[kms- 1 ] 


[kms" 1 ] 


[kms" 1 ] 


[kms- 1 ] 


[M ] 


[M Q ] 


[M ] 


[Myr] 

i 1+0.5 



66811 

190429 A 

15570 

14947 

210839 

163758 

192639 



40.0 
39.0 
38.0 
37.0 
36.0 
34.5 
33.5 



3.64 
3.62 
3.51 
3.52 
3.54 
3.41 
3.42 



5.91 
5.96 
5.94 
5.83 
5.80 
5.76 
5.68 



-5.70 
-5.68 
-5.66 
-5.85 
-5.85 
-5.80 
-5.92 



2300 
2300 
2200 
2300 
2100 
2100 
1900 



0.9 
1.0 
1.1 
1.3 
1.0 
1.1 
1.3 



0.05 
0.04 
0.05 
0.03 
0.05 
0.05 
0.05 



100 

30 
30 
30 
120 
30 
30 



210 
150 
97 
130 
210 
94 
90 



90 
57 
40 
36 
80 
34 
43 
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M' . refer to spectroscopic mass, initial mass and mass at the corresponding stellar age, respectively. Uncertainties for T cff are : 



^spec > 

1000 K, ± 0. 1 dex for log g and ±0.1 dex on log L/L Q (reflecting the uncertainty on the distance). An uncertainty of : 
estimated for the mass-loss rates, v™ is measured within ± 100 kms -1 . 
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Table 4. He, C, N, O and P abundances. 



Star 


y(He/H) 


e(C) 


e(N) 


e(O) 


e(P) 


e (E CNO) 


X(C) 


X(K) 


X(0) 


HD 16691 


0.15 


6.52±0.18 


9.00±0.17 


7.83±0.25 


5.05±0.24 


9.03 


1.00E-04 


8.63E-03 


6.67E-04 


HD 66811 


0.16 


7.60±0.25 


9.10±0.17 


8.13+0.30 


5.05±0.26 


9.16 


2.86E-04 


1.05E-02 


1.30E-03 


HD 190429 A 


0.15 


7.09±0.18 


8.92±0.16 


7.83±0.26 


5.26±0.20 


8.96 


9.10E-05 


7.27E-03 


6.69E-04 


HD 15570 


0.10 


7.52±0.16 


8.62±0.15 


8.30±0.24 


5.06±0.18 


8.81 


3.27E-04 


4.79E-03 


2.63E-03 


HD 14947 


0.12 


8.30±0.16 


8.78±0.12 


8.13±0.24 


5.25±0.15 


8.97 


1.66E-03 


5.00E-03 


1.44E-03 


HD 210839 


0.12 


8.22±0.21 


8.70±0.15 


8.48±0.14 


4.70±0.18 


8.99 


1.32E-03 


4.67E-03 


3.23E-03 


HD 163758 


0.15 


8.56±0.16 


8.78±0.16 


8.36±0.19 


5.02±0.21 


9.08 


2.72E-03 


5.19E-03 


2.25E-03 


HD 192639 


0.15 


8.17±0.17 


8.76±0.15 


8.61±0.21 


5.07±0.17 


9.05 


1.09E-03 


5.01E-03 


4.01E-03 


OriOBl (B stars) 




8.35±0.03 


7.82±0.07 


8.77±0.03 




8.92 


1.96E-03 


6.77E-04 


6.43E-03 


Sun 3D 


0.09 


8.39±0.05 


7.78±0.05 


8.66±0.05 


5.36±0.03 


8.88 


2.16E-03 


6.18E-04 


5.35E-03 



By convention log e(X) = 12 + log [X/H]. For comparison, we give the solar chemical composition determined by Asplund et al. (2005, Sun 
3D) and the chemical composition of B stars of the Ori OB 1 association determined by Simon-Diaz (2010) and Nieva & Simon-Diaz (201 1). 
Mass-fractions are also indicated for CNO. 



The luminosities were derived for the sample stars under the 
assumption that they are single. If they turn out to be binaries, the 
luminosity of each component would have to drop, thus translat- 
ing into lower initial masses and a somewhat different ages for 
the stars. As already stated (Sect. 2.1), seven of the eight stars 
in our sample are most likely single and the known binary is 
well separated, therefore we were able to use photometry of the 
primary for the luminosity calculation. 

When studying the most massive stars in the Arches clus- 
ter, Martins et al. (2008) found a significant overlap between the 
positions of the faintest WN7-9h stars and the most luminous, 
"extreme", O supergiants in their sample, which turned out to 
have spectral types 04-6If + . The positions of the stars of our 
sample in the H-R diagram overlap with the stars studied by 
Martins et al. (2008) of the same spectral class. We note that 
HD 190429A and HD 16691 show a late-type WN spectrum in 
the infrared K band (Conti et al. 1995) and could be regarded as 
"transitional" objects in terms of evolutionary properties, much 
like the extreme OIf + found by Martins et al. (2008), although 
they are systematically less luminous. When merged with those 
from Martins et al. (2008), our results confirm that the position 
of WN7-9h stars and O supergiants is distinct in the H-R dia- 
gram. 

Another insight into the evolutionary status of the stars is 
shown in Fig. 4, where the hydrogen mass fraction is plotted 
against the stellar luminosity. It is clear from Fig. 4 that all stars 
in our sample are on the early part of their evolution away from 
the zero-age main-sequence, even the most extreme of them 



( HD 190429A). Compared to the results from Martins et al. 
(2008) our stars appear to be more evolved (more hydrogen- 
depleted or equivalently helium-enriched) than their sample of 
supergiants, either normal or extreme. The reason for this differ- 
ence is unclear, because helium lines in the near-infrared are as 
sensitive to variations of helium content as lines in the optical 
range. A younger age of these cluster stars and possibly a higher 
metallicity could influence both their evolution and spectra. 

6.2. Masses 

The spectroscopic masses were derived after correcting the ef- 
fective gravity derived from the analysis by the effect of cen- 
trifugal forces caused by rotation (see Table 3), following the 
approach outlined in Repolust et al. (2004). Assuming a ran- 
dom distribution of orientations for the rotational axis, the true 
surface gravity g c is the sum of the effective gravity g e ff and a 
correction term. To a good approximation, the centrifugal accel- 
eration averaged over the stellar disk is given by the projected 
rotational velocity, whence g c = g e ff + (v sin i) 2 /R* 
For the stars in our study, the corrections are small, with 
Alogg < 0.04, which corresponds to a correction lower than 
10% for the actual mass. 

Evolutionary masses are obtained by bilinear interpolation 
between the tracks and isochrones in the H-R diagram. In Table 
3, we list masses at the corresponding stellar age and the origi- 
nal mass. Spectroscopic masses and evolutionary masses (at the 
star age) generally agree well. The difference is below 10% for 



12 



Bouret et aL: Properties of Galactic early-type O-supergiants 




4.7 4.6 4.5 4.4 

log(Teff) 



Fig. 3. H-R diagram showing the program stars (black squares) and the 
stars from Martins et al. (2008) for comparison. Blue triangles are "nor- 
mal" 04-6 supergiants, green squares are "extreme" 04-6If + stars, pink 
diamonds are WN7-9h objects. Evolutionary tracks and isochrones are 
from Meynet & Maeder (2003) with an initial rotation of 300 km s -1 . 



HD 16691, £ Puppis, HD 190429 A, HD 14947, HD 163758, and 
HD 192639, while it is around 20% for HD 15570 and A Cep. 
For the most part, the observed differences could result from un- 
certainties in log g. Assuming the stellar distance is known, the 
error in log g is the dominant source of error in the spectroscopic 
mass, hence in the ratio of the evolutionary mass to the spectro- 
scopic mass. These uncertainties are partially attributable to the 
difficulties associated with rectifying echelle spectra, especially 
in the vicinity of broad H lines. If accurate surface gravities are 
to be obtained, it is crucial that observing procedures are adopted 
to facilitate rectification, and to allow the accuracy of the recti- 
fication to be verified (e.g., measurement of multiple standards). 
Ideally, it would be useful to check the rectification against lower 
resolution (but still high-signal-to-noise) spectra. Obviously, dis- 
tance determinations to Galactic O stars are quite uncertain (at 
least for field stars) and we cannot discard just yet the contribu- 
tion of distance uncertainties as an important potential source of 
errors in spectroscopic masses. 

Likewise, continued progress in stellar evolution theory con- 
firms that stellar rotation is another fundamental physical pro- 
cess that affects the evolution of massive stars. The global agree- 
ment reached with spectroscopic masses furthermore suggests 
that, along with the introduction of NLTE line-blanketed stellar 
atmospheres, the assumptions adopted in the evolution models 
provide a realistic description of the actual conditions in young 
early-type massive objects. 

£ Puppis is very interesting in this context. It is likely seen 
near equator-on (Harries 2000, and refs. therein), which im- 
plies that its true rotation rate should be 210 kms~' or slightly 
higher. Applying the centrifugal correction to the measured sur- 
face gravity, we derive an actual (spectroscopic) mass of 56 M Q , 
while the initial mass, as derived from the star's location in the 
H-R diagram is 62 M Q (cf. Table 3). The evolutionary model 



with an initial mass of 60 M and an initial rotation rate of v sin i 
= 300 km s _1 that best represents the position of f Puppis shows 
that after 3.2 Myr, the stellar mass should be * 50 M G and the 
rotation rate m 108 kms -1 , more than 100 kms~' lower than the 
measured equatorial velocity of the star. Either the initial rota- 
tion velocity of £ Puppis was higher than 300 kms -1 or the star 
has been spun up during its evolution. Estimating the real value 
for the initial rotation rate is not possible because too many other 
quantities are involved in the evolution of the star's angular mo- 
mentum. For instance the mass-loss rates, as measured from this 
study, are a factor two to three times lower than the values from 
Vink et al. (2000), which were adopted in the evolution mod- 
els of Meynet & Maeder (2003). After 3.2 Myr a 60 M G (initial 
mass) star should have lost only 3 M G rather than 8 M Q (as- 
suming the other parameters would still evolve according to the 
higher mass-loss scenario). 




0.00 0.0028 0.0056 0.0084 0.0112 
log L/Lq 



Fig. 4. Hydrogen mass-fraction as a function of luminosity. Symbols 
are the same as in Fig 3. Full lines are for Geneva models with an initial 
rotational velocity of 300 km s -1 (Meynet & Maeder 2003). The dashed 
line indicates the initial value X(H)=0.705 



6.3. Surface abundances 
6.3.1. Helium and nitrogen 

To assess the chemical evolutionary status of the stars in our 
sample, we compared the nitrogen to helium mass-fraction ratio 
as a function of effective temperature to predictions from stellar 
evolution models. According to stellar evolution, nitrogen and 
helium are produced during the hydrogen burning phase but the 
relative increase of abundances is larger for nitrogen than for 
helium, since the latter is already abundant in the atmosphere 
while the former is not. The N/He ratio quickly evolves from the 
initial value to a value corresponding to the CNO equilibrium, 
and stars for which the equilibrium has not been reached will 
present intermediate values of N/He. 

Two of the 04 supergiants, HD 16691 and f Puppis, seem 
slightly offset in this diagram, and masses higher than inferred 
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from the H-R diagram would be needed to account for their po- 
sition. In both cases, this is a consequence of the strong nitrogen 
enrichment we derive rather than from too small a helium mass- 
fraction, which we find to be normal for this range of masses and 
ages. The case of HD 15570 is the exact opposite: a moderate ni- 
trogen enrichment and low helium abundance shifts the star up 
in the diagram. Although HD 190429 A is a binary, its properties 
seem well described by standard evolution models. The location 
of these two stars in Fig. 5 also indicates that they will soon reach 
the CNO equilibrium (possibly the case for HD 14947 as well). 

Similarly, the location of the later-type stars (cooler T e g) in 
this diagram is consistent, both qualitatively and quantitatively, 
with the predicted N/He fraction from stellar evolution models 
for objects in the mass range 50 - 40 M Q . These stars lie in a 
region of the plot where the N/He ratio is still rising, which in- 
dicates that the CNO equilibrium has not been attained. 
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Fig. 5. Ratio of nitrogen-to-helium mass-fraction as a function of effec- 
tive temperature, as predicted in Geneva models for stars with an initial 
rotation rate of 300 kms~' (Meynet & Maeder 2003), and as derived 
from this study (symbols). The red part of each line corresponds to the 
part of evolutionary tracks where X(H) > 0.4. 



6.3.2. Carbon versus nitrogen 

Together with the surface enrichment of nitrogen, the stars in our 
sample reveal carbon depletion, resulting in significant increases 
of the N/C ratios. This ratio is an excellent tracer of the chemi- 
cal evolution of massive stars, because it increases substantially 
when CNO cycle-processed material appears at the stellar sur- 
face. 

Figure 6 presents a comparison between the measured N/C 
ratios (mass-fractions) and the predictions of the Geneva stel- 
lar evolution models, as a function of stellar age. We can use 
this plot, together with the H-R diagram, to check whether a star 
is evolving according to the prediction of standard evolutionary 
tracks. We find two different behaviors - one for 04-type super- 
giants and the other for later-type (05 to 07.5) supergiants. The 
measured N/C ratio of this later group behaves qualitatively and 



quantitatively as predicted by stellar evolution models. These 
stars are less massive and older (around 4 Myr or older) than 
the first group of objects between 3 and 3.5 Myr. In the early- 
type, younger, group, N/C ratios are higher than predicted for 
the corresponding masses. ForHD 16691 and HD 190429 A, it is 
about ten times higher than predicted by any evolutionary model, 
which additionally confirms the advanced evolutionary status of 
these objects. 
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Fig. 6. Evolution of the ratio of nitrogen-to-carbon mass-fractions as a 
function of time in Geneva models (Meynet & Maeder 2003). The full 
(black) line is for a model with 60 M Q and an initial rotation velocity 
of 500 km s . Overplotted in red, green, blue and pink are the parts of 
the tracks corresponding to the main-sequence phase (more exactly the 
part of the hydrogen burning phase where X(H) > 0.4). 



Maeder et al. (2009) recently highlighted that N/C ratios de- 
pend on many variables including stellar age, mass, metallic- 
ity, rotational velocity, and multiplicity or magnetic fields and 
warned against an overinterpretation of measurements compared 
to theoretical predictions. However, taken at face value, these 
plots suggest that the global efficiency of mixing in stars more 
massive than 60 M Q is underestimated in stellar evolution mod- 
els for a given rotation rate. 

In Fig. 6 we also show the expected N/C ratio as a function 
of age for a 60 M Q model with an initial rotation velocity of 
500 kms -1 . The primary effect of increasing the initial rotation 
velocity is to maintain the N/C ratio at its maximum value up 
to 4.6 Myr (shift by 0.15 Myr compared to a model with 300 
kms" 1 ). This rotation increase has no effect in the region of the 
plot where the 04 supergiants are located, i.e. around 3 Myr. 
Hence, a spread in rotation rates appears unable to account for 
the higher N/C ratios we observe. 

Two stars of our sample, HD 66811 and HD 210839, have 
high rotation rates with vsin i > 210 kms 1 . HD 16691 has a 
lower vsin i (« 135 kms -1 ) but presents several spectral char- 
acteristics usually associated with a high rotation rate (N m and 
C m, He ii ,14686, Ho-, see Sect. 5.3 and Appendix A). HD 668 1 1 
and HD 16691 have the strongest nitrogen enrichment of the 
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sample and HD 16691 the lowest carbon abundance (hence the 
highest N/C ratio). Both stars support the idea that rotational 
mixing is responsible for the surface composition changes. On 
the other hand, HD 210839, which also has a high rotation rate, 
goes against this conclusion because it has the lowest nitrogen 
enrichment and marginal carbon depletion. The other stars have 
intermediate v sin i and no correlation emerges with the N/C ra- 
tio. This contradictory result is reminiscent of the findings of 
Hunter et al. (2007), though we stress that our sample is too 
small for a statistical interpretation. 

Let us now compare how the carbon and nitrogen mass- 
fractions relate to each other in a diagram where X(C) is plot- 
ted as a function of X(N). It is expected that in the early phases 
of stellar evolution, X(C) decreases, while X(N) is built up at 
the expense of carbon. Fig. 7 shows that, as expected, stars with 
higher nitrogen content show stronger carbon depletion. For the 
measured stellar masses, and for stars still evolving in the hy- 
drogen burning phase, the observed range of the nitrogen mass- 
fractions is qualitatively compatible with the predictions of stel- 
lar evolution models. However, quantitative differences are ob- 
served. Six out of eight stars cluster in a region where the carbon 
depletion is too strong for the measured nitrogen enrichment (or 
alternatively, nitrogen enrichment is too low for the measured 
carbon depletion). Once again, we observe a dichotomy in the 
behavior of the 04 supergiants and the 05 - 07.5 supergiants. 
The former group of stars, with masses « 60 M Q , shows the 
strongest carbon depletion, which is so extreme in two cases 
that it would not be reached in any model at any evolution- 
ary phase. In this context it is interesting that these stars show 
WN spectral features in their near-infrared spectra (Conti et al. 
1995). Only HD 14947 (and HD 163758 to a lesser extent, given 
its very peculiar abundances) agree qualitatively and quantita- 
tively with the 40 - 60 M Q models that are expected to best de- 
scribe the properties of this sample. This agreement probably 
means that both stars are actually very close to CNO equilib- 
rium. Furthermore, the position of HD 163758 in Fig. 7 is con- 
sistent with its carbon-rich status as quoted from observations by 
Gomez & Niemela (1987b). 

6.3.3. CNO 

In massive stars, the total mass fraction of CNO remains roughly 
constant during the CNO cycles. In the early phases of evolution, 
the CN-equilibrium is achieved first while the oxygen content 
remains about constant, and somewhat later enters the burning 
cycle. The full CNO-equilibrium is achieved only after a signifi- 
cant amount of hydrogen has been burnt in the core. Changes in 
the oxygen surface abundance are only expected in later stages 
in the more massive and faster rotating models. 

In Fig. 8, we compare our N/C versus N/O ratios to those pre- 
dicted by stellar evolution models for different masses and X(H) 
> 0.4 (Meynet & Maeder 2003). This plot indicates whether the 
variations of N/C and N/O follow the trend expected for material 
processed through the CNO cycle. The amplitudes of the depar- 
ture from the initial abundance ratios are notoriously sensitive to 
differences in rotational velocities, masses, and ages. 

For the 04 supergiants, N/C ratios are too high for the mea- 
sured N/O ratios, compared to predictions by the evolution mod- 
els. For three stars, their location in Fig. 8 indicates that their N/C 
and N/O ratios are consistent with more advanced evolutionary 
status (WNs) than actually derived. The fourth star, HD 15570, 
exhibits different properties: although its N/C and N/O ratios are 
slightly higher than those predicted by evolutionary models for 
60 M e , they are nevertheless compatible with the its evolution- 
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Fig. 7. Carbon mass-fraction as function of nitrogen mass-fraction, as 
predicted by Geneva models with an initial rotation velocity of 300 
kms~' (Meynet & Maeder 2003). Solid line is for 85 M Q , dashed for 
60 M Q , and dot-dashed for 40 M e . Symbols are for values derived from 
this study. 



ary status. Within this group of stars, the observed behavior is 
qualitatively compatible with the expected evolution: the more 
massive stars are farther along the evolutionary tracks than the 
less massive ones, assuming that they are in the same evolution- 
ary phase/status (as suggested by the H-R diagram, cf. Fig. 3). 

The opposite behavior is observed for later- and mid-type 
supergiants, whose measured N/C ratios are lower than predicted 
from evolution models for stars in the 40 - 60 M Q range (for the 
measured N/O ratios). 

The case of HD 163758 is noteworthy because it presents 
unusual abundance patterns (see Sect. 5.7), including a super- 
solar carbon abundance and a very high nitrogen content while 
being only mildly oxygen-depleted. Such patterns likely result 
from a non-standard (non solar-like) original mixture of CNO 
elements, and the relative enhancement of the nitrogen mass- 
fraction depends on the initial carbon abundance available for 
conversion during the CNO-cycle. 

Figure 7 indicates that the nitrogen enrichment in HD 14947 
is consistent with the carbon depletion predicted for a 50 M G star. 
Moreover, Fig. 6 indicates that the measured N/C ratio is consis- 
tent with the age of the star. From these two results, we conclude 
that the measured N/C ratio is compatible with the predictions of 
stellar evolution models. Therefore, the position of HD 14947 in 
Figure 8 suggests that N/O is too high compared to evolutionary 
predictions, or in other words, that the measured oxygen abun- 
dance is too low (lower than predicted). If the oxygen depletion 
is not the result of abnormal initial abundances, it could be ex- 
plained if deeper layers of the star are exposed to the surface, for 
example by deeper mixing than predicted by the rotating mod- 
els. Alternatively it could be explained by stripping of the star 
by winds or a binary companion; or by the accretion of oxygen- 
depleted layers from a former companion (Langer et al. 2008; de 
Mink etal. 2011). 
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We note that the total of the CNO abundances of the sample 
stars (see Table 4) shows that with a mean combined abundance 
e (2CNO) = 9.02, the initial CNO mixture for several objects 
is potentially different from the assumed standard cosmic abun- 
dances adopted for the computation of the stellar evolution mod- 
els by Meynet & Maeder (2003). Because there is evidence for 
a significant spread in the metallicity for solar-type stars in the 
thin disk (<x[Fe/H] = 0.191 dex, Fuhrmann 2011), there might 
be a spread in the actual Fe abundances in O-type stars as well. 
Although solar values were assumed in our study, they might not 
be relevant for all our sample stars. 

Although striking, the scatter in N/C and N/O ratios for 
this sample is real, as indicated by their very similar spectra 
(hence physical parameters) but different CNO line intensities. 
For the first group of stars (HD 16691, HD 1 90429 A, HD 668 1 1 , 
HD 15570), the amplitude of the deviation from the models is 
low enough that differences in the initial abundances could help 
reconcile the outputs from stellar evolution models and our mea- 
surements. The slope of the theoretical tracks is very sensitive to 
the adopted initial mixture (see, e.g., Fig. 3 of Przybilla et al. 
2010) and a span of N/O and N/C ratios is expected. Including 
magnetic fields in stellar evolution models could also change the 
speed and amplitude at which abundance patterns are built up at 
the surface. A recent study by Meynet et al. (2011) confirmed 
that the inclusion of magnetic braking yields very different re- 
sults for the surface abundances, which depend on the assumed 
rotation law inside the star. For stars with differential rotation, 
mixing is faster and stronger, but for stars with solid-body ro- 
tation, mixing is inhibited because the stars spin down rapidly. 
In this case, surface abundance ratios are lower than in models 
without magnetic braking. However, Meynet et al. (2011) ob- 
tained results for a 10 M Q star and dedicated simulations would 
be necessary to confirm that more massive stars with magnetic 
fields can still be rotating relatively normally and produce sig- 
nificant changes in surface abundance ratios after a few Myr. 

For the second group of (later) stars, an obvious solution ap- 
peals to models with higher rotation velocities, which exhibit 
more mixing. Again, if magnetic braking plays a role, mix- 
ing would be modified (see above). Alternatively, (close) bi- 
narity would affect the evolution and likely the surface abun- 
dances of the stars (from mixing and possibly mass transfer). 
Binarity would also imply different luminosities, hence evolu- 
tionary masses. In any case, comparison with tailored models 
for magnetic braking (and/or lower masses) would be required 
for more quantitative conclusions. 

The interplay between rotational mixing and mass-loss as 
currently implemented in the models might not be efficient 
enough to account for the observed the CNO patterns. We note 
in this context that we derive mass-loss rates up to three times 
lower than those used in the calculations of evolution models 
(see Sect. 6.4). The actual metal content in the Geneva evolu- 
tionary models is also higher than what is measured in the so- 
lar neighborhood (0.02 versus 0.014) by Przybilla et al. (2008) 
(2008), which best corresponds to the global metallicity of the 
stars of the present sample. Both these facts undermine the direct 
comparison of our measurements with chemical patterns from 
stellar evolution models for an assumed Z/Z Q = 0.02. 

The overall properties of these objects, including their loca- 
tion in the H-R diagram and hydrogen mass-fraction (see Fig. 
3 and Fig. 4), indicate that they are observed at different stages 
of their chemical evolution, even though they are of very similar 
age. The underlying drivers of the observed differences are un- 
clear and could be manyfold - from enhanced mixing related to 
binary evolution in an earlier phase prior to the supergiant stage 



to evolution with rotation and a magnetic field (Meynet et al. 
2011). Recent results by Sana & Evans (2011) suggest that the 
binary fraction among massive stars is at least as high as 50%. 
Although all but one star in the present sample are believed to 
be single, we cannot exclude the possibility that some of them 
might have been members of a binary system in the past; this 
scenario is especially appealing for runaway stars like f Puppis 
(or HD 16691). On the other hand, the observed properties of 
the only massive binary in our sample, HD 190429A, are qual- 
itatively similar to those of the other stars of the same spectral 
type, and no clear influence of the secondary in this system is 
found. Nevertheless, since binarity is expected to have signif- 
icant consequences for stellar evolution, tailored evolutionary 
models would be needed to compare to our measurements and 
to draw firmer conclusions. 
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Fig. 8. N/O vs. N/C abundance ratio (by mass) measured for the sample 
stars and compared to the ratios predicted in Geneva models (Meynet 
& Maeder 2003). Solid line is for 120 M , dashed red line is for 85 M , 
dot-dashed green line is for 60 M Q , and dotted blue line is for 40 M Q . 
The black dot-dot-dot-dashed line is for a model with 60 M Q and v sin i 
= 500 km s _I . Only the part of evolutionary tracks corresponding to the 
evolutionary phase where X(H) > 0.4 are plotted. 



6.4. Wind properties 

6.4.1. Mass-loss rates and filling factors 

Clumping as considered in this paper is often referred to as 
micro-clumping because it assumes that the clumps are optically 
thin. The filling factors we derive cluster around / M = 0.05 with 
a fairly limited dispersion +0.02. This corresponds to mass-loss 
rates reduced by factors of four to six compared to the homoge- 
neous values that would be derived from fitting Ha. More inter- 
estingly, for all stars our mass-loss rates are within a factor of 
three of the predictions by Vink et al. (2000). Fig. 9 displays the 
ratio of predicted to measured mass-loss rates as a function of 
effective temperature. The ratios seem to increase with T e g (the 
same trend is observed with luminosity). Theoretical mass-loss 
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Fig. 9. Ratio of theoretical calculations by Vink et al. (2000) to mea- 
sured mass-loss rates as function of effective temperature. Theoretical 
mass-loss rates are computed for the best-fit parameters for each star. 
The uncertainties in the ratios are estimated using error bars for mea- 
sured M only. 



A promising alternative comes from the work by Lucy (2007, 
2010a,b), who showed that turbulence in the photosphere re- 
duces the wind driving force at the sonic point, thus strongly 
reducing the mass outflows, especially when compared to the 
calculations of Vink et al. (2000). Calculations were made in 
the microturbulence framework, because the nature and cause of 
the actual turbulence are unknown. It is expected that microtur- 
bulence should depend on the effective temperature, hence the 
reduction in mass-loss rates (relative to Vink et al. 2000) should 
be higher at higher T e g, which might reduce the trend observed 
on Fig. 9. 

Most likely, a proper description of stellar winds from mas- 
sive stars should include optically thick clumps together with 
non-void interclump medium, and a non-monotonic velocity 
field (Sundqvist et al. 2011). This is especially important in ac- 
counting for velocity-porosity effects, which Sundqvist and col- 
laborators conclude is especially relevant for resonance lines. 
Using a complex mixture of ID radiation-hydrodynamics mod- 
els to create 2D and 3D stochastic wind models, they computed 
synthetic line profiles for A Cep, based on a Monte-Carlo radia- 
tive transfer code. They were able to reproduce the observed pro- 
files of P v resonance lines and Ha and to a lesser extent of He n 
A4686. They derived a mass-loss rate of 1.5 x 10~ 6 M Q yr _1 , 
which is much higher than the mass-loss that they would de- 
rive assuming optically thin clumps only (4. x 10~ 7 M G yr _1 ), 
although this last value is obtained from the fit to the P v reso- 
nance lines only. 



rates by Vink et al. (2000) decrease for decreasing effective tem- 
perature for fundamental stellar parameters typical of early-mid 
O-type supergiants. Therefore, Fig. 9 shows that the theoretical 
mass-loss rates increase more rapidly with T e fi (alternatively L) 
than measured in the present sample. 

As r e ff increases, the maximum of the flux distribution grad- 
ually shifts to shorter wavelengths, where more lines are avail- 
able to contribute to the driving, thus increasing the mass-loss 
rates. The difference apparent in Fig. 9 suggests that more 
driving is present at higher T e $ in the theoretical calculations. 
However, we used a very large number of species and ions for 
our calculations of hydrodynamics in the models, at least as 
many as were included in Vink et al. (2000). Since mass-loss 
rates are mostly set by ions of the iron group, it is unlikely that 
we missed important lines that would instead be accounted for 
in the theoretical calculations by Vink et al. (2000). 

Small-scale clumping significantly alters the ionization bal- 
ance in the wind (Hillier 1991; Bouret et al. 2005). Depending 
on the temperature regime, dominant ions for different species 
might change, which could lead to a direct effect on the radia- 
tive driving. A recent study by Muijres et al. (201 1) argued that 
the impact of clumping on the wind driving is an intricate func- 
tion of the stellar temperature, the density inside the clumps and 
the physical size of the clumps. In the specific case of optically 
thin clumps, increased recombinations lead to an increase in the 
mass-loss rate, as a consequence of an increased line force. The 
amplitude of this effect is stronger as the clumping factor in- 
creases. This would be somewhat compensated for if porosity 
effects were included in their calculation, as a consequence of 
a decreased interaction between photons and the gas. All things 
considered, empirical mass-loss rates such as those we derived 
here cannot be reconciled as yet to those resulting from theoret- 
ical consideration, even including clumping corrections. 



The mass-loss rate we derive for A Cep is very close to the 
one quoted by Sundqvist et al. (201 1). Our fit to the wind lines is 
equally good when adopting lower-than-solar phosphorus abun- 
dance. Interestingly, both works fail to produce a double-peaked 
emission profile at He n /14686. It turns out that this specific 
shape, (as well as the absorption trough seen in Ha) is quite 
sensitive to the actual treatment of the rotation law in the wind. 
Sundqvist et al. (2011) had to adopt a fairly large amount of 
macro-clumping to reproduce the central absorption in Ha. The 
extent to which a more proper treatment of the fast rotation of 
this star would change this conclusion is unclear. We will present 
models of A Cep and £ Puppis including a realistic treatment of 
rotation in a forthcoming paper, but we argue at this step that 
this ingredient is as important as clumping (either optically thin 
or thick) to explain the observed shapes of Ha and He n /14686 
and that neglecting it leads to overemphasize the effect of (large- 
scale) clumping. At the same time the volume-filling factor ap- 
proach slightly overestimates the strength of Ha (Sundqvist et al. 
201 1). All in all, A Cep might not be as representative of the in- 
fluence of wind structures on emerging spectra as initially ex- 
pected. 

Clumping was found to start close to the photosphere in all 
but three stars our of sample, two of which are fast rotators. For 
the third one, namely HD 16691, we already mentioned that sev- 
eral emission profiles are characteristic of fast rotation, although 
not reflected by the vsin i. We realize that the need to start 
clumping farther out in the wind is mostly driven by the need 
to fit the absorption trough seen in Ha in the case of f Puppis 
and A Cep. Such a trough is not as clearly seen in HD 16691 (as 
M is much higher) and the fit to N iv /11718 was the driver for the 
adopted v cl . In any case, an accurate treatment of rotation is re- 
quired for these objects to better investigate the actual properties 
of clumping. 
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6.4.2. Optically thin clumps and the P v problem 

We showed in previous studies (Bouret et al. 2003; Hillier et al. 
2003; Bouret et al. 2005) that the volume filling factor associ- 
ated with wind clumping can be derived from the simultaneous 
fit to P V-U1118-1128, O v ,11376, N iv ,11718, He n ,14686, and 
Ha. The abundances for helium, nitrogen and oxygen were fixed 
from the analysis of photospheric lines; hence we had no other 
leverage but clumping and mass-loss rates to fit the wind lines 
listed above. The case of phosphorus is somewhat different be- 
cause the FUV resonance doublet is the only available feature. 
Since the abundance of phosphorus is not expected to change 
over the lifetime of an individual star, a solar abundance should 
be appropriate. However, we had to use lower-than-solar phos- 
phorus abundances, in addition to clumping, to improve the fit to 
P v resonance doublet. This result is reminiscent of our findings 
(see Hillier et al. 2003; Bouret et al. 2005) and although it might 
hint at a selective depletion of phosphorus in O supergiants with 
respect to the other metal abundances, it is also possibly related 
to the assumptions adopted to describe clumping in the models. 
The recent introduction of velocity porosity in radiative trans- 
fer calculations by Sundqvist et al. (2011) showed that the P v 
AAl 118-1 128 lines of A Cep can be fitted with a solar abun- 
dance of phosphorus. We note in this context that these lines in 
HD 163758 do not reflect the ratio of oscillator strengths of the 
two lines, thus indicating that these lines are saturated. Prinja & 
Massa (2010) concluded from their study of Si iv ,1,11393-1403 
on a sample of B supergiants that velocity porosity could explain 
this behavior. HD 163758 is therefore a very good candidate to 
investigate a more realistic treatment of wind inhomogeneities 
for O supergiants. 

The primary consequence of being forced to reduce the phos- 
phorus abundance in CMFGEN models is to change the radial 
optical depth of the resonance lines, hence their strength. Since 
the optical depth of the wind in one line is related to the ion- 
ization fraction of the corresponding ion (r ra j oc Mq^A^ with 
Ae and q{ the abundance of element E and its ionization frac- 
tion for stage i), the need to decrease the phosphorus abundance 
compared to solar could be related to the incorrect determina- 
tion of the wind ionization (i.e., over- production of the P +4 ion). 
To investigate the problem of phosphorus ionization in more de- 
tails, we derived an empirical Mq(P +4 ) by profile fitting using 
the Sobolev with exact integration (SEI) method (Lamers et al. 
1987). We refer to Massa et al. (2003) for a detailed presenta- 
tion of the method. We used wind quantities such as Voo and /3 
from the spectroscopic study (Table 3). Figure 10 compares the 
Mq(P +4 ) derived with this approach to those from the modeling 
with CMFGEN. Overall, the agreement between the quantities 
from both techniques is very satisfactory when a reduced phos- 
phorus abundance is used. In other words, in this case only, very 
similar optical depths are derived, while when P = P is used, 
the optical depths from CMFGEN are too high. Although the 
global scaling with P/P Q < 1 improves the agreement, a clear 
trend for too high Mq(P +4 ) in CMFGEN is observed at high ve- 
locities (typically above v/Voo> 0.7) for all but one star. The fit 
(Appendix A) shows that the synthetic CMFGEN profiles are in- 
deed too strong in absorption near the blue-edge of the absorp- 
tion components of the P Cygni profiles of P v lines 5 . This di- 
vergence between Mq(P +4 ) from CMFGEN and SEI is stronger 

5 f Puppis shows by far the largest disagreement for empirical and 
theoretical Mq(P +4 ), which also translates to the poorest fit to the ob- 
servations within the sample. The clear peak in the SEI Mq(P +4 ) distri- 
bution between 0.6 < v/v«,< 0.7 might indicate the presence of a DAC 
in the wind of the star at the time of the observation 



for stars hotter than 36,000 K. It indicates that the ionization 
fraction of q(P +4 ) is predicted too high in CMFGEN in the outer 
part of the winds. The ionization fraction of P +4 in our models 
shows that P v is a dominant ion below 0.5v M for all stars, con- 
sistent with results in Fullerton et al. (2006) and Rrticka & Kubat 
(2009). In stars with 7/ eff > 37,000 K, P v remains the dominant 
phosphorus ion throughout the wind. 

Sundqvist et al. (2011) suggested that X-rays have a signif- 
icant influence on the ionization in the outer part of the winds. 
We performed several tests with different X-ray fluxes for these 
stars. We did not find a significant influence of X-ray radiation 
on the ionization of phosphorus, unless very high X-ray fluxes 
were used, which would alter the ionization of other species as 
well. Such very high X-ray fluxes provides a "tail" of XUV radi- 
ation that is strong enough to alter the P v ionization fraction but 
the corresponding (emitted) logLx/Lboi are not consistent with 
the now well-defined relation (Sana et al. 2007). We note that 
Rrticka & Kubat (2009) already concluded that X-ray radiation 
alone cannot change the ionization fraction of P v significantly 
enough to avoid M reduction when clumping is taken into ac- 
count. 
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Fig. 10. Mq(P +4 ) as a function of normalized velocity derived from the 
modeling with CMFGEN (black) or from SEI (red). 



On the other hand, Waldron & Cassinelli (2010) concluded 
that strong wind XUV-EUV emission line radiation could alter 
the ionization fraction of P v, while remaining compatible with 
observational constraints. Nevertheless, other important ques- 
tions such as photoionizations from excited states, influence of 
other species and optical depth effects in the He n continuum that 
might challenge their conclusion are ignored. 

Our fits to the P v AAl 1 18-1 128 resonance lines are very sat- 
isfactory, given the relative simplicity of the adopted treatment 
of clumping. Again, the adopted phosphorus abundance is lower 
than solar but we seem to need the same P/P = 0.5 ± 0.2 ratio 
for all stars, which we can adopt as a fiducial value for sub- 
sequent modeling of other early-type supergiants. We believe 
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that wind parameters as derived from the present study, by fit- 
ting several wind profiles over a range of species and ioniza- 
tion stages provide a realistic and robust description of the stel- 
lar winds. Nevertheless, more consistent descriptions of wind 
clumping will need to be considered in future modeling. The 
velocity porosity presented by Sundqvist et al. (201 1) is promis- 
ing, although the velocity structure of the hydrodynamic models 
do not exhibit the porosity in velocity space that is required to 
produce the observed P-Cygni profile shapes, and reconcile Ha 
and UV mass-loss rates. Other important questions will need to 
be addressed, including a self-consistent calculation of the ion- 
ization structure, the extension to optically thick winds, and a 
realistic description (as yet lacking) of the structure of stellar 
winds. 



6.5. Line force consistency check 

Since CMFGEN computes the full radiation field, we can accu- 
rately compute the radiative force. Consequently, we can check 
whether the derived radiative force is consistent with that needed 
to drive the wind. For convenience we write the momentum 
equation as 



dv 1 dP„ 

grad = V— + - — +g 

dr p dr 



(4) 



where P g is the gas pressure, v is the velocity, g = GM/ r 2 , and 
grad is the total radiative force, which is evaluated using 



grad 



Ancpr 2 



(5) 



whereby is the flux mean opacity. 

If the mass-loss rate, velocity law, and gravity in our models 
are consistent, the left-hand-side (LHS) of Eq. 4 should be the 
same as the right-hand-side (RHS). In Fig. 11 we compare the 
LHS and RHS for a model similar to that found for £ Puppis: 
r eff =41,000K, M= 1.7 x lO^Moyr 1 , and f x = 0.05 (which 
can be compared to our best-fit model of T eB =40,000 K, M- 
2.0 x 10~ 6 M G yr _1 , and /«, = 0.05). As is readily apparent, our 
line force is too strong in most of the wind, indicating that the 
velocity law is incorrect, or more likely, the mass-loss rate is too 
low. The same conclusion is reached if we use the parameters 
found for the best-fit £ Puppis model. 

The simplest means of obtaining better consistency between 
the LHS and RHS would be to increase the mass-loss rate, and 
the simplest way of maintaining "consistency" with the spectral 
fit would be to simultaneously increase the volume filling fac- 
tor /. As expected, tests show that a moderate increase in the 
mass-loss rate to 2.7 x 10~ 6 M Q yr~' and faa =0.1 yields a much 
better agreement. As discussed previously, a slightly higher / 
value might give consistent spectral fits if the effects of velocity 
porosity were taken into account. 

Because we can perform this consistency check, we could, 
in principle, use the wind dynamics to derive an additional con- 
straint on M. However, there are several difficulties involved. 
First, it is very difficult to obtain agreement between the LHS 
and RHS around the sonic point (v ~ 20 km s -1 ) - this is not 
surprising, since the critical point lies near the sonic point. We 
can only obtain a moderate agreement if we adopt a low micro- 
turbulent velocity of 5 km s _1 , a value which is inconsistent with 
that required to fit the Fe lines in the UV, and with that gener- 
ally inferred for O supergiants. Understanding the nature of this 
microturbulence is of crucial importance for understanding radi- 
ation driven winds. Moreover, we found it easier to achieve a fit 



when /3 < 1 in this region, a value more in line with the predic- 
tions of radiation-driven wind theory. Second, there is increasing 
evidence that winds are clumped and dynamic - not steady-state 
as assumed by Eq. 4. Therefore, for full consistency, we need 
to run hydrodynamic simulations, and since £ Puppis is a rapid 
rotator, this may necessitate 3D simulations. Third, there is the 
problem of velocity porosity, which will tend to reduce the line 
force for a given mass-loss rate and clumping. Fourth, we need to 
simultaneously reduce the abundances, because mass-loss rates 
are predicted to scale as Z 085 (Vink et al. 2001). A systematic 
error in the abundances by a factor of 1.5 (which is not unrealis- 
tic) will cause a ~ 40% error in M. Consequently, while we can 
conclude that the spectroscopically derived mass-loss rates, and 
the mass-loss rate inferred from the line force, are consistent to 
within a factor of 2 for £ Puppis, it is almost impossible to obtain 
better consistency at the present time. It is also worth noting that 
it is necessary to include Ne, S, Ar, CI etc. in the model calcula- 
tions. While they do not have a major influence on the observed 
spectrum, they do significantly enhance the radiative line force 
in the wind. 

In general, the derived spectroscopic mass-loss rates, and 
those inferred from the dynamics for O supergiant agree rea- 
sonably well. The same is not true for dwarf O stars, where 
the derived CMFGEN line force indicates a much higher mass- 
loss rate than that derived spectroscopically from the weak UV 
P Cygni profiles (e.g., Martins et al. 2012). 



7. Conclusions 

We have modeled the FUV-UV and optical spectra of a sample of 
eight O-type supergiants to investigate their surface abundances 
and wind properties. The results that we obtained can be sum- 
marized as follows: 

- Supergiants with spectral types between 04 to 07.5 order 
along a well-defined evolutionary sequence, in terms of ages 
and masses, from younger and more massive to older stars 
with lower initial masses. 04 supergiants cluster around the 
3 Myr isochrone and are more massive than 60 M Q , while 
from 05 to 07.5, stars have masses in the range 50 - 40 M Q 
and are 4±0.3 Myr old. No mass-discrepancy is observed 
since spectroscopic and evolutionary masses agree within the 
uncertainties (and in any case within 20%). Despite the high 
signal-to-noise ratios that can be achieved, the use of echelle 
spectra compromises the rectification of spectra, especially 
in the vicinity of the broad H lines that are the primary di- 
agnostics of log g. Improving log g determinations is crucial 
for determining stellar masses by spectroscopic methods. 

- All stars exhibit significant enrichment of nitrogen, as ex- 
pected from their evolutionary status. Carbon is depleted 
with respect to an initial abundance assumed to be solar-like 
for all but one star, namely HD 163758. This star turned out 
to be carbon-rich (twice solar at least) as expected from its 
spectral classification. If depleted by the amount expected 
on the basis of its evolutionary status, its initial carbon abun- 
dance was likely very high. All other stars have either a ni- 
trogen enrichment that is too low for the measured carbon 
depletion, or alternatively too strong a carbon depletion for 
the measured nitrogen enrichment. 

- The observed carbon and nitrogen mass-fractions are com- 
patible with those expected from the models, for the mea- 
sured stellar masses. On the other hand, the N/C ratios as 
a function of age are inconsistent with the theoretical pre- 
dictions for the four earliest (04 spectral type) stars of the 
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Fig. 11. Illustration comparing the radiative line force with other terms in the momentum equation. The model parameters are similar to those of 
£ Puppis; r c ff=41,000 K and M Q = 1.7 x 1CT 6 M Q yr~' . The right-hand plot illustrates the sensitivity of the line force to the adopted microturbulent 
velocity around the sonic point. Below 10km s -1 hydrostatic equilibrium is valid. 



sample. The efficiency of rotational mixing as a function of 
age is questioned for these stars. Other mechanisms might be 
needed to explain the observed patterns. A solution could be 
that the observed metallicity gradient within the galactic disk 
should be taken into account rather than assuming standard 
solar values for metals. Alternatively, a modified efficiency 
of mixing may be needed in stellar evolution models. 

- Mass-loss rates derived using clumped models are lower than 
theoretical predictions by factors that range from 1.5 to 2.5. 
This difference is actually very reasonable given the known 
uncertainty on both estimates. The corresponding filling- 
factors associated with small-scale clumping are 0.05+0.02. 
Clumping is found to start close to the photosphere for all but 
three stars, two of which are fast rotators, namely A Cep and 
£ Puppis. Although its vsin i is smaller, HD 16691 behaves 
like these two fast rotators; several optical lines indicate that 
this star could indeed be a fast rotator as well, seen at a higher 
inclination. 

- The ratios of theoretical to observed mass-loss rates increase 
with luminosity (or r en -), which suggests that too much driv- 
ing might be present in the models for the hotter and more 
massive stars. We found no convincing explanation for this 
behavior at this point. The need for a sub-solar phospho- 
rus abundance to fit the P v resonance line was investigated. 
Overall, empirical ionization fractions derived within the 
SEI framework agree well with the theoretical fractions from 
CMFGEN. Some remaining discrepancies are observed, es- 
pecially in the outer regions of the winds (higher radial ve- 
locities). P v is predicted to be overabundant by CMFGEN 
in these circumstances. In any case, sub-solar phosphorus 
abundance are needed to artificially decrease the radial op- 
tical depth of the P v lines. Although this problem is directly 
related to the micro-clumping framework we work with, the 
needed abundance reduction is fairly constant and could be 
adopted as default for future work. 

To reach a more complete, more quantitative, understand- 
ing of the observed trends of the photospheric and wind proper- 
ties for Galactic supergiants, we will extend our analysis to later 
spectral types (down to the end of the O-type stars sequence). In 
forthcoming papers, we will also explore the different technical 
problems we met during the present analysis, such as the forma- 
tion of excited lines of CNO elements that are fundamental for 



an abundance analysis; and the extent to which incorporating ro- 
tation improves fits to key diagnostic features in the spectra of 
massive stars. 
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Appendix A: Best fits 



In this appendix we present our best-fit models to the optical and 
UV spectra for the eight stars in our sample. The wavelength 
range between 1200 and 1225 A was not used in the spectral 
analysis since it suffers from a strong interstellar Lyman absorp- 
tion. Nevertheless, we take this interstellar absorption into ac- 
count by adding the effects of a representative hydrogen column 
density to our synthetic spectra in our final plots. 
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Fig. A.l. Best-fit model for HD 190429A (red line) compared to FUSE and WE spectra (black line) 






V = 




— v — ^_ 



(D 4440 




i u I 





AIN 

i 


1 






NV— CI II 

1 II 1 _ 











1 





Wavelength (A) 

Fig. A.2. Best-fit model for HD 190429A (red line) compared to the ELODIE spectrum (black line) 
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Fig. A.3. Best-fit model for HD 668 1 1 (red line) compared to Copernicus and IUE spectra (black line) 
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Fig. A.4. Best-fit model for HD 6681 1 (red line) compared to the FEROS spectrum (black line) 
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Fig. A.5. Best-fit model for HD 16691 (red line) compared to FUSE and IUE spectra (black line) 
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Fig. A.6. Best-fit model for HD 16691 (red line) compared to the ELOD1E spectrum (black line) 
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Fig. A.7. Best-fit model for HD 15570 (red line) compared to FUSE and IUE spectra (black line) 
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Fig. A.8. Best-fit model for HD 15570 (red line) compared to the ELOD1E spectrum (black line) 
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Fig. A.9. Best-fit model for HD 14947 (red line) compared to FUSE and IUE spectra (black line) 
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Fig. A.10. Best-fit model for HD 14947 (red line) compared to the ELODIE spectrum (black line) 
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Fig. A.ll. Best-fit model for HD 210839 (red line) compared to FUSE and 1UE spectra (black line) 
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Fig. A.12. Best-fit model for HD 210839 (red line) compared to the NARVAL spectrum (black line) 
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Fig. A.13. Best-fit model for HD 163758 (red line) compared to FUSE and IUE spectra (black line) 
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Fig. A.14. Best-fit model for HD 163758 (red line) compared to the UVES spectrum (black line) 
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Fig. A.15. Best-fit model for HD 192639 (red line) compared to FUSE and 1UE spectra (black line) 
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Fig. A.16. Best-fit model for HD 192639 (red line) compared to the ELODIE spectrum (black line) 



